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Black-hole (BH) binaries with single-BH masses m = 5-20M©, moving on quasicircular orbits, 
are among the most promising sources for first-generation ground-based gravitational-wave (GW) 
detectors. Until now, the development of data-analysis techniques to detect GWs from these sources 
has been focused mostly on nonspinning BHs. The data-analysis problem for the spinning case is 
complicated by the necessity to model the precession-induced modulations of the GW signal, and by 
the large number of parameters needed to characterize the system, including the initial directions of 
the spins, and the position and orientation of the binary with respect to the GW detector. In this 
paper we consider binaries of maximally spinning BHs, and we work in the adiabatic-inspiral regime 
to build families of modulated detection templates that (i) are functions of very few physical and 
phenomenological parameters, (ii) model remarkably well the dynamical and precessional effects 
on the GW signal, with fitting factors on average > 0.97, but (iii) might require increasing the 
detection thresholds, offsetting at least partially the gains in the fitting factors. Our detection- 
template families are quite promising also for the case of neutron-star-black-hole binaries, with 
fitting factors on average « 0.93. For these binaries we also suggest (but do not test) a further 
template family, which would produce essentially exact waveforms written directly in terms of the 
physical spin parameters. 

PACS numbers: 04.30.Db, x04.25.Nx, 04.80.Nn, 95.55.Ym 


I. INTRODUCTION 

A world-wide network of laser-interferometer gravitational-wave (GW) detectors, recently built 0, has by now 
begun operation. Inspiraling binaries of compact objects, such as black holes (BHs) and neutron stars (NSs) are 
among the most promising astrophysical sources for these detectors. The GWs from the inspirals are expected to 
enter the frequency band of good detector sensitivity during the last few seconds or minutes of evolution of the 
binaries; GW scientists plan to track the phase of the signals very accurately, and to enhance the signal-to-noise 
ratio by integrating the signals coherently over their duration in the detector band. This is achieved by filtering the 
detector output with a bank of templates that represent our best theoretical predictions for the signals. 

Until now, the development of data-analysis techniques has been focused mostly on binaries containing NSs (whose 
spins are negligible for data-detection purposes) and nonspinning BHs 0 . Nonspinning, high-mass BHs pose a delicate 
problem: the breakdown of the post-Newtonian (PN) expansion in the last stages of the inspiral makes it hard to 
prepare reliable templates for the detection of binary BHs (BBHs) of relatively high total mass (say, 10-40Mq) 
with LIGO-VIRGO interferometers. Various resummation techniques, such as Fade approximants Q and Effective 
One-Body (EOB) techniques Hi have been developed to extend the validity of PN formalism |^. Damour, Iyer, 
and Sathyaprakash 0 compared the templates generated by different PN treatments, and found that they can be 
very different. In a companion paper to the present one henceforth BCVI], we investigated this issue for the GW 
signals emitted by comparable-mass BBHs with a total mass M = lO-dOM©. In BCVI we proposed a few examples of 
detection template families (DTFs), built either as time series or directly in the frequency domain, which try to address 
the failure of the PN expansion. The philosophy behind DTFs is to replace a family of signals that correspond to a 
specific mathematical model of the binary with families that can cover a broader range of plausible signals. Because 
the direct correspondence with the mathematical model is lost, DTFs are appropriate for the purpose of first detecting 
GW signals, but do not give direct estimates of physical parameters, such as the masses of the binary constituents. 
[Within the EOB framework, see also the recent paper by Damour, Iyer, Jaranowski, and Sathyaprakash where 
the authors extend 3PN EOB templates with seven flexibility parameters and then show that the unextended 3PN 
templates already span the ranges of the flexibility parameters consistent with plausible 4PN effects.] 

Very little is known about the statistical distribution of spins for the BHs in binaries: the spins could very well 
be large. Apostolatos, Cutler, Sussman, and Thorne Hll henceforth ACST] have shown that when this is the 
case, the evolution of the GW phase and amplitude during the inspiral will be significantly affected by spin-induced 
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modulations and irregularities. In a BBH, these effects can become dramatic if the two spins are large and they 
are not exactly aligned or antialigned with the orbital angular momentum. If this happens, there is a considerable 
chance that the analysis of interferometer data, carried out without taking spin effects into accounts, could miss the 
signals from these spinning BBHs altogether. The gravitational waveforms from binaries of spinning compact objects 
depend on many parameters: the masses and spins of the objects, the angles that describe the relative orientations of 
detector and binary, and the direction of propagation of GWs to the detector. In practice it is impossible, due to the 
extremely high computational cost, to filter the signals with a template bank parametrized by all of these parameters. 
One strategy is that of providing effective templates that depend on fewer parameters, but that have still reason ably 
high overlaps with the expected physical signals. An interesting suggestion, built on the results obtained in Ref. |lfl| . 
came from Apostolatos |lll . who introduced a modulational sinusoidal term in the frequency-domain phase of the 
templates to capture the effects of precession. However, while Apostolatos’ family reduces the number of parameters 
considerably, its computational requirements are still very hig h. Moreover, using an approximated analytical model of 
NS-BH waveforms, Grandclemeiit, Kalogera and Vecchio showed that this family fails to capture those waveforms 
satisfactorily (see however Ref. for a hierarchical scheme that can improve the fit by adding “spikes?? in the 
template phasing). 

In this paper, complementary to BGVI, we study the data analysis of GWs from binaries with spinning BHs; for 
simplicity, we restrict our analysis to the adiabatic limit, where the two compact objects in the binary (either two 
BHs, or a NS and a BH) follow an adiabatic sequence of spherical orbits driven by radiation reaction (RR). The 
denomination of spherical orbits comes from the fact that the orbital plane is not fixed in space, but precesses, so 
the orbits trace a complicated path on a (slowly shrinking) spherical surface. We neglect the problems caused by 
the failure of PN expansion in these binaries (note that the conservative part of the EOB framework ^ has already 
been extended to the spinning case by Damour 0, providing a tool to move beyond the adiabatic approximation; 
we plan to add radiation-reaction effects to this model, and to study the consequences on GW emission and detection 
elsewhere). Here, we carry out a detailed study of PN processional dynamics and of GW generation in processing 
binaries in the adiabatic limit, and we use the resulting insights to build a new class of modulated effective templates 
where modulational effects are introduced in both the frequency-domain amplitude and frequency-domain phase of 
the templates. The mathematical structure of our templates suggests a way to search automatically over several of the 
parameters (in strict analogy to the automatic search over initial template phase in the data analysis of nonspinning 
binaries), reducing computational costs significantly. We argue that our families should capture very well the expected 
physical signals. 

We note here a shift in perspective from BCVl. In this paper, we use the PN equations for the two-body dynamics 
of spinning compact objects to build a fiducial model (our target model) that represents our best knowledge of 
the expected physical signals. Because we cannot use the target model directly for data analysis (it has too many 
parameters), we build effective template families with fewer parameters. These families are then compared with 
the target model for a variety of binary parameters, to gauge their ability to match the physical signals (their 
effectualness i). On the other hand, in BGVI we employed several variants of the PN equations (with diverging 
behaviors in the late phase of inspiral) to identify a range of plausible physical signals; we then built our DTFs so that 
they would match all of the PN target models satisfactorily. This said, we shall still refer to the template families 
developed in the present paper as DTFs. We direct the reader to BGVI for a simple introduction to matched-filtering 
techniques and their use in GW data analysis (developed in the literature by various authors Hill) , and for an 
explanation of some of the notation used in this paper. 

This paper is organized as follows. In Sec. m we define the target model, and we explain the conventions used to 
represent the generation and propagation of GWs. In Sec. IHIl we study the two-body dynamics of spinning compact 
objects, looking for the features that are especially relevant to the data-analysis problem. Using this insight, in 
Sec. lIVI we formulate our DTFs, and we also describe two families of standard stationary-phase-approximation (SPA) 
templates, to be used as a comparison when evaluating the performance of the DTFs. In Sec. El we discuss the 
overlap and false-alarm statistics of our DTFs. In Sec. EH we evaluate the performance of our DTFs for BBHs and 
NS-BH binaries, and we briefly discuss a more advanced (and very promising) template family for NS-BH systems. 
In Sec. Em we summarize our conclusions. 

Throughout this paper we adopt the noise spectral density for LIGO-I given by Eq. (28) of BGVI. The projected 
VIRGO noise curve is quite different (deeper at low frequencies, with a displaced peak-sensitivity frequency). So our 
results for high-mass binaries cannot be applied naively to VIRGO. We plan to repeat this study for VIRGO in the 
near future. 
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II. DEFINITION OF THE TARGET MODEL 

In this section we define the target model used in this paper as a fiducial representation of the GW signals expected 
from processing binaries of spinning compact objects. We restrict our analysis to the adiabatic regime where the 
inspiral of the compact objects can be represented as a sequence of quasicircular orbits. At any point along the 
inspiral, a binary of total mass M = mi + m 2 and symmetric mass ratio rj = miTO 2 /M^ is completely described by 
the orbital angular frequency lo, the orbital phase 'I', the direction Ljv oc r x v of the orbital angular momentum, and 
the two spins Si = yimiSi and S 2 = X 2 m 2 S 2 , where Si ,2 are unit vectors and 0 < xi ,2 < 1- Throughout this paper 
we shall use carets to denote unit vectors, and we shall adopt geometrical units. 

In Sec. Ill Al we write the PN equations that govern the adiabatic evolution of the binary and the precession of Ljv 
and of Si, 2 - All the target waveforms used to test the effectualness 0 of our DTFs are obtained by integrating these 
equations numerically. The validity of the adiabatic approximation is discussed in App. ^ In Sec. urm we discuss 
our criterion for stopping the numerical integration of the evolution equations at the point where the adiabatic 
approximation ceases to be valid. In Sec. lira building on Refs. i we describe a formalism for computing 

the response of a ground-based detector to the GWs generated by a spinning binary; the response is not just a 
function of the trajectory of the binary, but also of the relative direction and orientation of binary and detector. The 
formalism describes also how the precession of the binary modulates the detector response. Last, in Sec. Ill d 1 we give 
a classification of all the parameters that enter the expression for the detector response, distinguishing those that 
specify the evolution of the binary itself from those that describe the relative direction and orientation of binary and 
detector. 


A. Equations for an adiabatic sequence of processing spherical orbits 


The path of the binary across the sequence of quasicir cular orbits is described by the adiabatic evolut ion of th e 
orbital angular frequency to up to 3.5PN order |lA l2ri7| with spin effects included up to 2PN order |2lll7ll(i| . 
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where tk = 0.577... is Euler’s constant, and where 9 is an arbitrary parameter that enters the GW flux at 3PN 
order and that could not be fixed in the regularization scheme used by the authors of Ref. j^. Note that in 
Eq. (Q we set the static parameter Wg = 0 (Note for v4 of this paper on gr-qc: Eq. (1) is now revised as per 

Ref. the parameter 9 has been determined to be 1039/4620 li^.l T he precession equations for the two spins are 
(see, for instance, Eqs. (4.17b,c) of Ref. or Eqs. (11b,c) of Ref. [loj| 'l 
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where we have replaced r and |Liv| by their leading-order Newtonian expressions in to, 

|Liv| =nr^uj = 77M^/3^-i/3. 
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This approximation is appropriate because the next spin-precession term is higher than the leading order, 

while next terms in the expressions of r and |Liv| are higher. 
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The precession of the orbital plane (defined by the normal vector L^v) can be computed as follows. From Eqs. (4.7) 
and (4.11) of Ref. 0 we see that the total angular momentum J and its rate of change Jrr (due to RR) depend on 
w, Ln and Si ^2 (schematically) as (S = Si + S 2 ) 


J = L + S = ?7M2(Ma;)-i/3LN \l + - ry (Mw)^/^ Ses+S , 


(5) 
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is known as effective spin M- Note that both terms in the L brace of Eq. © originate from orbital angular momentum 
(the second term comes from the spin-orbit coupling). Taking the time derivative of we obtain 
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where to get the last term on the right-hand side we have used uj = Comparing Eqs. ||HJ| and projecting 

out only the direction perpendicular to Ln, and keeping only the terms up to the leading and next-to-leading orders, 
we get 
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Thus, we now have the set of four equations ©-© and @ for the four variables w. Si, S 2 , and Ln- We follow 
Ref. |l(i|. Eq. (4.15), in defining the accumulated orbital phase 'k as 
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( 10 ) 


This phase describes the position of the two compact objects along the instantaneous circular orbits of the adiabatic 
sequence; the phase of the GW waveforms, as detected by a ground-based detectors, differs from this by processional 
effects, as explained below in Sec. ITTcl 


B. Endpoint of evolution 


The orbital energy of the two-body s ystem at 2PN and 3PN orders, expressed as a function of w, and assuming the 
static parameter Wg = 0 |2^ , reads Il7l2ll2fl 
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( 12 ) 


In the context of our adiabatic approximation, it is natural to stop the integration of Eqs. 0-0 and at the point 
(the Minimum Energy Circular Orbit, or MECO) where the energy i?„pN reaches a minimum. 
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after this point the adiabatic approximation breaks down [^. (The MECO is discussed by Blanchet for non¬ 
spinning binaries under the name ICO, for Innermost Circular Orbit.) However, if we find that tL> = 0 (which implies 
certainly that the adiabatic approximation has become invalid) before the MECO is reached, we stop the evolution 
there. In BCVl we noticed that for nonspinning binaries this behavior occurs for the 2.5PN evolutions, but not at 
2PN, 3PN and 3.5PN orders. 

Throughout this paper, we shall call the instantaneous frequency of GWs at the endpoint of evolution the ending 
frequency, which, up to a correction that arises from processional effects, is twice the instantaneous orbital frequency 
defined in this section. It so happens (see BCVl) that a knowledge of the ending frequency is important to cut off 
the candidate detection templates at the point where we know too little about the physical signals to model them 
further. In Sec. IIII 5l we study the dependence of the ending frequency on the spins of the binary. 


C. Gravitational waveforms 


As we have seen, the trajectory of the inspiraling binary is obtained by integrating Eqs. 0-0 and 0 for the time 
evolution of w(t), Si(t), S 2 (t) and hpfft). To determine the corresponding gravitational waveforms, we need to choose 
a specific coordinate system. We follow the convention proposed by Finn and Chernoff la henceforth EC], and also 
adopted by Kidder EC employ a fixed (source) coordinate system with unit vectors {ef, e®, ef} (see Fig. 0. 
For a circular orbit, the leading-order mass-quadrupole waveform is (throughout this paper, we use geometrical units) 




(14) 


where D is the distance between the source and the Earth, and where Q)? is proportional to the second time derivative 
of the mass-quadrupole moment of the binary. 


= 2 [\^ n^] , (15) 

with n* and A* the unit vectors along the separation vector of the binary r and along the corresponding relative 
velocity v. These unit vectors are related to the adiabatic evolution of the dynamical variables by 

n = ef cos $5 -I- ef sin , A = — ef sin -|- ef cos <i>s ; (16) 

the vectors ef 2 form an orthonormal basis for the instantaneous orbital plane, and in the EC convention they are 
given by 
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FIG. 2: 


Detector and radiation frames in the FC convention M- 


The vector ef points in the direction of the ascending node of the orbit on the {x, y) plane. The quantity is the 
orbital phase with respect to the ascending node; its evolution is given by 

$ 5 =w —dcosi, (18) 

where i and a are the spherical coordinates of Ln in the source frame, as shown in Fig.^] Using Eqs. d and llTCll . 
we can write Eq. d as 

= - 2 ([e^]'^ cos2$s + sin 2 ^ 5 ), (19) 

where the polarization tensors and e'^ are given by 

e;^ = ef (g) ef — ef 0 ef , ef = ef 0 ef + ef 0 ef . (20) 

Eor a detector lying in the direction N = ef cos0 + ef sin0, it is expedient to express GW propagation in the 

radiation coordinate system with unit vectors {ef^,ef,ef-} (see our Fig. ^ together with, for instance, Eq. (4.22) of 
Ref. given by 

COS0 —ef sin0, (21) 

= ef, (22) 

ef = ef sin 0 + ef cos 0 = N . (23) 

In writing Eqs. we used the fact that for a generic binary-detector configuration, the entire system consisting 

of the binary and the detector can be always rotated along the z axis, in such a way that the detector will lie in the 
{x, z) plane. Later in this paper (in Sec. nvii we shall find it convenient to conserve the explicit dependence of our 
formulas on the azimuthal angle tp that specifies the direction of the detector. 

In the transverse-traceless (TT) gauge, the metric perturbations are 

h^^ = /i+T+ + /ixTx, (24) 

where 
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The response of a ground-based, interferometric detector (such as LIGO or VIRGO) to the GWs is M 
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Binary 


GW propagation 

Detector orientation 

M, rj. Si, S2 

9si , 6s2 , <^Si - <t>S2 

5 0Sl 0S2 

0, ifi 

9, (j>, if 

Basic 

Local 


Directional 



TABLE I: Classification of binary, GW-propagation, and detector parameters. 


where and Fx are the antenna patterns, given by 
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with ex,y the unit vectors along the orthogonal interferometer arms. For the geometric configuration shown in Fig.|21 
with detector orientation parametrized by the angles 9, (f> and ip, we have 


F+ = -(1 + cos^ 0) cos2(/) cos2^ — COS0 sin20 sin27/>, 

Fx = -(1 + cos^ 0) cos2(/) sin2'0 + COS0 sin2(/) cos2^ . 

Inserting Eqs. ini), and (123 into Eq. 123, we get the final result 0. 

^resp = Cq cos2$5 + Sq sin2$5 , 
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Sq = [S+F+ + SxFx], 
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(37) 


D. Binary and detector parameters 


We shall refer to the total mass M, to the mass ratio rj = mi TO 2 /M^, and to the magnitudes of the two BH (or NS) 
spins Si and S 2 as the basic parameters of the binary. Once these are set, we complete the specification of a binary 
configuration by giving the initial orbital phase and the components of the orbital and spin angular momenta in the 
source frame, for a given initial frequency. In our convention, the initial orbital angular momentum is determined by 
the angles (0 Ln = L'/'Ln = Q^), Eis shown in Fig.|3 [In the rest of this paper, we shall use the notation (0 Ln,<('Ln) to 
emphasize that these angles define the direction of the orbital angular momentum.] The directions of the spins are 
specified by the angles {Os^, ipSi) and (^Sj, ^^ 52 ), defined with respect to an orthonormal basis aligned with Ln, 


_ Ln X ef _ ^ ^ 

ei = -— , 02 = Ln X ei, 03 = Ln , 

|Ln X 0*1 

shown in Fig. 01 We then have 

51 = 01 sin^Sj cos 4 >Si +62 sin^Sj sinc^s^ +03 cos^Sj , 

52 = 61 sin 9 s 2 cos ((>52 + 02 sin 0^2 sin + 03 cos 632 ■ 


(38) 


(39) 

(40) 










FIG. 3 : Specification of the initial Newtonian orbital angular momentum in the source frame {ex,ey,ez}. 



FIG. 4 : Specification of the initial directions of the spins with respect to the FG orthonormal basis {61,02,03} [Eq. lEHl]. 


Among the six angles (0 Ln j ^Ln ) j (^Si j 4>Si ) j and , (j)s 2 ) > only three are intrinsically relevant to the evolution of the 
binary: ^ 5 ^, Os^ and (j)Si — 4’S2- We shall refer to them as local parameters. The other three independent parameters, 
which are relevant to the computation of the waveform, describe the rigid rotation of the binary as a whole in space, 
and we shall refer to them as directional parameters. In fact, there are five more directional parameters: 0 and 
(/? specify the direction to the detector in the source frame, and 0 , cj), and tp specify the orientation of the detector 
with respect to the radiation frame. All these parameters have already been introduced in the previous section. Our 
classification of the 15 binary and detector parameters is summarized in Tab.^ 









9 


III. ANALYSIS OF PRECESSIONAL DYNAMICS 

In a seminal paper M , ACST investigated in detail the evolution of binaries of spinning compact objects, focusing 
on orbital precession and on its influence on the gravitational waveforms. In this section, we build on their analysis 
to discuss several aspects of quasicircular precessional dynamics that are especially important to the formulation of 
a reliable DTFs for these systems. Note also that Wex has derived analytic solutions for quasielliptical solutions 
to the 2PN conservative dynamics, including spin-orbit effects. 

We complement ACST’s analytical arguments with the empirical evidence obtained by studying the orbits generated 
by the numerical integration of Eqs. and We select the following typical binaries: BBHs with masses 

{nil + 1 ^ 2 ) given by (20 -I- 1O)M0, (15 -I- 15)M0, (20 -I- 5)Mq, (10 -I- 1O)M0, (7 -I- 5)Mq-, and NS-BH binaries with 
masses mi = IOMq (BH) and m 2 = IAMq (NS). The BHs are always chosen to be maximally rotating {S = m?)^ 
while the NSs are assumed to be nonspinning. There are neither astrophysical data nor theoretical results which 
suggest that maximal spins are preferred. However, in this paper we decide to investigate the most pessimistic (in 
terms of precessional effects) scenario. The initial GW frequency is chosen at 30 Hz for binaries with total mass larger 
than 20 Mq, and 40 Hz for all the other cases. For each set of masses, we consider 1000 (or, when the numerical 
study is very computationally expensive, only 200) orbital evolutions obtained with random initial orientations of the 
orbital and spin angular momenta. (These initial configurations are taken from the pseudorandom sequence specified 
in Sec. ME and used in Sec. iwn to evaluate the effectualness Q of our DTF in matching the target signals.) 

In Sec. imsi we introduce the ACST results, and in particular the distinction between simple and transitional 
precession. In Sec. ME we study the dependence of the GW ending frequency (defined in Sec. Ill Bll on the initial 
values of spins and on their evolution, and we link this dependence to the conservation of certain functions of the 
spins through evolution. As mentioned above, a knowledge of the ending frequencies of our target model is important 
to decide what extension each of the detection templates should have in the frequency domain. In Sec. inm we 
examine the value of the binding energy and of the total angular momentum at the end of evolution, and we estimate 
the amount of GWs that must be emitted during plunge, merger and ringdown to reduce the spin of the final BH 
to the maximal value. In Sec. Em we discuss, largely on the basis on numerical evidence, the effects of spin on 
the accumulated orbital phase 4* [defined by Eq. (llOll ]: we argue that these effects are mainly nonmodulational, and 
that, for data-analysis purposes, they can be treated in the same way as the standard PN corrections to the orbits 
of nonspinning binaries. It follows that the precession of the orbital angular momentum is the primary source of 
modulations in the signal (as already emphasized by ACST for particular classes of binaries). In Sec. IIII hi we show, 
again on the basis of numerical evidence, that transitional precession has little relevance to the data-analysis problem 
under consideration. In Sec. Em we discuss the power-law approximations introduced by ACST to describe the 
precession of the orbital plane as a function of frequency in particular binaries, and we show that they are appropriate 
in general for the larger class of binaries under consideration in this paper. These approximations are a basic building 
block of the effective template families developed by Apostolatos and, indeed, of our generalized and improved 
families. 


A. The ACST analysis 

In their paper M on precessing binaries of compact objects as GW sources, ACST chose to work at the leading 
order in both the orbital phasing and the precessional effects to highlight the main features of dynamical evolution. Eor 
orbital evolution, they retained only the first term on the right-hand side of Eq. O: as a consequence, the precession 
of the orbital plane is the only source of GW modulation considered in the analysis. [The resulting accumulated orbital 
phase 4 /, given by Eq. is known as Newtonian Chirp.] Eor the precession of the orbital angular momentum and 
of the spins, ACST retained only the first terms (the spin-orbit terms) in Eqs. ®, ® and ®. On the basis of 
these approximations, and in the context of binaries with either mi ~ m 2 (and spin-spin terms neglected) or S 2 ~ 0, 
ACST classified the possible evolutions of spinning binaries into two categories: simple precession and transitional 
precession. 

The vast majority of evolutions is characterized by simple precession, where the direction of the total angular 
momentum remains roughly constant, and where both the orbital angular momentum and the total spin S = Si + S 2 
process around that direction. ACST provided a simple analytic solution for the evolutions in this class. They also 
showed that the orbital precession angle, expressed as a function of the orbital frequency, follows approximately a 
power law (see Eq. (45) of Ref. [13 )■ 

Transitional precession happens when, at some point during the evolution, the orbital angular momentum and the 
total spin become antialigned and have roughly the same magnitude, so the total angular momentum is almost zero. 


J = L-kS « 0. 


(41) 
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FIG. 5: Binary ending frequencies (gray dots) as functions of the initial value of Kbs/k^^, for 1000 initial spin configurations 
of M = (15 + 15)M0 BBHs (in the left panel), and M = (10 + 1.4)M0 NS-BH binaries (in the right panel), at 2PN and 3.5PN 
orders. The solid lines plot the SO-only predictions for the ending frequencies. 


When this condition is satisfied, the total angular momentum is liable to sudden and repeated changes of direction. 
The evolutions in this class cannot be easily treated analytically, but they occur only for a small portion of the possible 
initial conditions. 

In this paper, we shall refer to the special cases investigated by ACST (with either toi « TO 2 , or S 2 ~ 0) as 
ACST configurations. NS-BH binaries and BBHs with mi 3> m 2 are astrophysically relevant cases among ACST 
configurations, because for both we can set S 2 ~ 0. The ACST formalism can also describe well BBHs with equal 
masses, but where spin-spin effects are negligible. 


B. Conservation laws and GW ending frequencies 


For the ACST configurations, both the total spin and its projection on the orbital angular momentum are constants 
of the motion: 


LN(i)-S(t) 


J ACST 


[S\t)] 


ACST 


const; 
const. 


(42) 

(43) 


For generic non-ACST configurations (as discussed, for instance, by Damour M), the effective spin Ses [Eq. Q] 
can, to some extent, replace the total spin in these conservation laws. 

From Eqs. and ©, we see also that if we ignore the spin-spin effects in the precession equations, then the 

projection 


_ Ln • Seff 
= M2 

of the effective spin onto the Newtonian orbital angular momentum is a constant of motion, 


(44) 


[Keff(t)]so = const (45) 

(where the subscript “SO” stands for the inclusion of spin-orbit effects only); on the other hand, neither S^(t) nor 
Sgjf(t) are conserved. 

The conservation of Kes has important consequences for the endpoints of evolution, defined in Sec. lHBi bv Eq. 11311 
for the MECO. In the nonspinning case, as discussed in BCVl, if the dynamics was known at all PN orders, then the 
MECO would agree with the Innermost Stable Spherical Orbit (ISCO), defined as the orbit beyond which circular 
orbits become dynamically unstable. When only spin-orbit (henceforth, SO) effects are included, the conservation 
of Keff preserves this correspondence between MECO and ISCO, because the leading-order SO term in the energy 
is proportional to Keg: in fact, the frequency of the MECO has a precise functional dependence on Keff [see Eqs. 
((TT |i -fTl ^ ]. 

When spin-spin (henceforth, SS) couplings are also included, KeS is no longer conserved, and the MECOs (and 
therefore the ending frequencies) of binaries with the same initial KeS become smeared around their SO-only values. 
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FIG. 6: Binary ending frequencies (gray dots) as functions of the initial value of for 1000 initial spin configurations 

of M = (20 + 5)Mq BBHs, at 2.5PN order. The solid lines plot the SO-only predictions for the ending frequencies. 


which are functions only of Keff. In addition to this smearing, the SS contribution to the energy introduces also a 
bias. In the end, however, the SS correction is not very important for the ending frequencies, as we can see in the 
following examples. In the left panel of Fig. 0 we plot the ending frequency at 2PN and 3.5PN orders versus 
the initial value of Keff for BBHs with M = (15 + 15)Mq (in gray dots), as compared to the SO-only predictions (in 
solid lines). The smearing of the ending frequencies is relatively mild, and so is the systematic deviation from the 
SO-only predictions. We have checked that this behavior characterizes all the mass configurations enumerated just 
before Sec. mrei at both 2PN and 3.5PN orders. In the right panel of Fig. we plot the ending frequencies for 
NS-BH binaries [with M = (10-1- 1.4)Mq]. The ending frequencies follow exactly the expected functional dependence 
on Keff. 

The mildness of these deviations can be understood (in part) by looking at the variation of Keff during the evolution. 
For example, for the (15 -I- 15)M0 BBHs shown in Fig. 0 the maximum deviation of Keff from being a constant 
(measured as maxdev(Keff) = [max(Keff) — min(Keff)]/2) is 0.036, to be compared with the maximum kinematically 
allowed deviation, 0.875; for (2O-|-5)M0 BBHs at 2PN order, maxdev(Keff) = 0.033, to be compared with the maximum 
kinematically allowed deviation 0.92. 

As we can infer from Fig.|3 the ending frequencies depend also on the PN order, and the difference between 2PN and 
3.5PN orders is more striking for NS-BH binaries than for BBHs. This trend is present also in the nonspinning case 
(see BCVl): for nonspinning {xi = X 2 = 0) equal-mass BBHs, we have = 0.137M“^ and = 0.129M“^. 

To give a few numbers, for a (10-1- 1O)M0 BBH, we have /gwSpn = ^43 Hz, and /qw^pn = 416 Hz; for a (15-1- 15)M0 
BBH, /gwSpn = 295 Hz, and /qwSpn = 277 Hz; on the other hand, for a (10 -I- 1.4) Mq NS-BH binary, we have 
/gw^pn = "^34 Hz, and = 559 Hz. For the second and third binaries, these values can be read off from the 

solid lines of Fig. 0 by setting Keff = 0 (no spins). 

Finally, in Fig. Owe show the ending frequencies for (20 -I- 5) Mq BBHs, when Eq. (0) (which rules the evolution of 
the orbital phase) is evaluated at 2.5PN order. In this case, if Keff > 0.5, then 6j goes to zero before the MECO can 
be reached. The resulting ending frequencies deviate considerably from SO-only predictions. As already discussed in 
BCVl, oj goes to zero because at 2.5PN order the gravitational flux goes to zero for high orbital velocities; since this 
very nonphysical behavior happens systematically, we then choose to exclude the 2.5PN order from our analysis. 


C. Energy radiated during inspiral and (estimated) total angular-momentum emitted after inspiral 

It is interesting to evaluate how much energy can be radiated in GWs before the final plunge, especially for binaries 
whose inspiral end in the LIGO-VIRGO frequency band. In the left panel of Fig. [ 7 | for M = (15 -I- 15)M0 BBHs, we 
plot the ratio between the 2PN (nonrelativistic) energy, given by Eq. (Hill and evaluated at the endpoint of evolution 
(as defined in Sec. irm . and the total mass-energy initially available, M. Depending on the initial relative orientation 
between the spins and the orbital angular momentum (as expressed by the initial Keff/K^’'), the energy that can be 
released in GWs during the inspiral ranges between ^ 1.5 and 3.5% of M. More energy can be emitted when the initial 
spins are aligned with the orbital angular momentum. We find similar results for all the other BBHs investigated, 
and similar results were also obtained by Damour in the EOB framework (see Eig. 1 of Ref. [l^). 
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FIG. 7: For 1000 (15 -I- 15)M0 BBHs with different initial spin configurations, in the left panel we plot the ratio between the 
(nonrelativistic) 2PN energy [Eq. lllll l at the ending frequency and the mass-energy initially available M, versus the initial 
in thfi right panel we plot the ratio between the total angular momentum J at 2PN order and the square of the 
(relativistic) 2PN energy [Eq. lllll l at the ending frequency, versus the initial Ksb/k^^. 


It is also interesting to estimate how much total angular momentum can be radiated during the coalescence phases 
that follow the inspiral (plunge and merger), especially when those phases fall in the LIGO-VIRGO band. In general, 
we have 


Jrad — J S 


BH ; 


(46) 


where Jrad is the angular momentum radiated during plunge-merger, J is the total angular momentum of the binary 
at the end of the inspiral, and Sbh is the spin of the final black hole. A lower limit on the angular momentum radiated 
in these phases can be obtained using the fact that the magnitude of the final spin can be at most Mgjj (where Mbh 
is the mass of the final black hole). To derive this lower limit we follow Flanagan and Hughes |^, and we write, 
using Eq. (ESI) . 

IJradI > |j| - ISbhI > |j| - Min > |J| - El, , (47) 


where Erei = M + E is the relativistic energy of the binary at the end of inspiral; in deriving Eq. lITTI) we used the 
relation |Sbh| < -^bh — -^rei- straightforward to see from Eq. I 47 II that this lower limit is nontrivial (that is, 
greater than zero) only when jJj > 

In the right panel of Fig.|7| for M = (15-|-15)M0 BBHs, we plot [Jj/if^gB where the angular momentum is evaluated 
at 2PN order |2lllT^ . 


J/M^ = 




^(81 - 577 ,+ 7 , 2 ) 


1 / Si 


7 , \M^ IVP 


-3 L 


M2 


M2 


M2 


S 

JP 


(48) 


We see that J is generally less than one, except when KeS > 0.4 (which happens for 13% of all the initial spin 
configurations); the maximum value of |J|/if2gj is 1.13. (For a similar plot obtained in the EOB framework see Fig. 
2 of Ref. El .) Such large values of KeS imply large ending frequencies [for the (15 + 15)M0 BBHs shown, larger 
than 400 Hz], which do not lie in the LIGO-VIRGO band of good interferometer sensitivity, unless the BBHs have 
higher masses; then all the frequencies are scaled down. In any case, for Kett = 1 (spins and orbital momenta initially 
aligned), in the high-mass binaries investigated, Eq. I 47 II suggests the lower limit 


IJradI >0.130.1 m2, (49) 

to be compared with the value 0.4M2 obtained by Flanagan and Hughes [ll using BH spins aligned with the orbital 
angular momentum (estimated to be ~ 0.9M2). 

A (trivial) upper limit for J^ad is obtained by setting Sbh = 0: 

IP'^'^I < |J|. 


(50) 
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Maximum modulational correction 


(2O + 1 O)M 0 (15 + 15)Mq (2O + 5 )M 0 (10 + 10)Me (7 + 5 )M 0 (10 + 1 . 4 )M 0 [NS-BH] 


(Ad '^®®)200 

0.0247 

0.0214 

0.0450 

0.0402 

0.0828 

0.1228 

^^90% (200) 

0.0460 

0.0411 

0.0676 

0.0787 

0.1504 

0.1884 

max 2 oo Atf®® 

0.0680 

0.0523 

0.1227 

0.1186 

0.2196 

0.1895 


TABLE II: Maximum modulational effects in the accumulated orbital phase We give the average over the 200 samples, the 
90% quantile of the distribution, and the maximum value for the diagnostic defined in Eq. 151-il . 


For different values of Keffj the upper limit for our (15 + 15)M0 binary is ^ 0.5 - I.IM^. However, in order for the 
inspiral to end within the LIGO-VIRGO band of good interferometer sensitivity (which requires a MECO frequency 
lower than 400Hz), we need KeS < 0.4, which corresponds to an upper limits ^ 0.5 - 0.7M^. 

To put this section into context, we point out that most reliable PN estimates for the energy and the angular 
momentum radiated after the MEGO can be achieved only with models that include information about the plunge 
phase, such as the model that can be built on Damour’s spinning-EOB equations HI- 


D. Spin-orbit and spin-spin effects on the accumulated orbital phase 

While for nonspinning binaries the accumulated orbital phase [defined by Eq. 11 OH ] coincides with (half) the GW 
phase at the detector, for spinning binaries the two phases differ by precessional effects; in the EG convention, these 
are found in part in the relation 


4 ) 5 . = 4 / — d cos i, 


(51) 


where is the orbital phase with respect to the ascending node of the orbit, which appears in Eq. m for the 
detector response to GW; and in part in the explicit time dependence of the coefficients Cq and Sq on a and i [see 
Eqs. I32ll - I37|l ]. In this section, we are going to argue that the evolution of the accumulated orbital phase is very 
similar in spinning and nonspinning binaries; and that, as a consequence, the effect of spins on detector response 
through the accumulated orbital phase can be reproduced using nonspinning-binary templates, such as those studied 
in BCVl [see also Eqs. Of course, precessional effects do enter the detector response through the other 

dependences mentioned above, and these cannot be neglected when building templates to detect physical signals. 

Both the spin-orbit and spin-spin couplings can affect the accumulated orbital phase dl through the 1.5PN and 
2PN terms in Eq. (Q. However, as we shall discuss in this section, this effect is largely nonmodulational. For each 
binary configuration, we introduce three different functions of time: (a) the accumulated orbital phase obtained 
by solving the full set of Eqs. 0-0 and ©, including the SO and SS couplings; (b) the accumulated orbital phase 
obtained by using the initial orbital angular momentum and spins at all times in the SO and SS couplings; and 
(c) the accumulated orbital phase for a nonspinning binary, obtained by dropping the SO and SS couplings 

altogether. 

In general, and are quite different for the same set of binary masses. However, the difference — 

yjynospin jg ^ strongly oscillating function (that is, it does not show any modulation), and it can be reduced 
considerably by modifying the I.5PN and 2PN coefficients in the phasing equation for the nonspinning binary. It is 
then reasonable to assume that such a nonmodulational effect could be captured by the nonspinning DTFs constructed 
in BCVl. Moreover, the difference between and is due to the nonconservation of the SO and SS terms that 
appear in Eq. d) for ilj. These terms have relatively high PN orders, so we expect that they will be small. 

Thus, we expect that 41^“^^ can be well described by a nonmodulational phasing of the kind 


d- 


nonmod 


(/) -Co + Cif + , 


(52) 


which looks rather like the frequency-domain phasings employed in the DTFs of BCVl. [Here C 2 and C 3 can be 
seen as actual (intrinsic) template parameters, whereas Cq and Ci represent, respectively, the initial phase and the 
time of arrival of the GW signal, both of which are extrinsic parameters in the sense discussed in BCVl.] To verify 
this hypothesis, we first evaluate in the frequency range 50 Hz-250 Hz (which is appropriate for first-generation 
ground-based GW detectors), using Eqs. 0-@ and @ at 2PN order, for all the BBH and NS-BH configurations 
considered earlier [(5 -I- 1) masses x 200 angles]. We then (least-square) fit with functions of the form (I52II . A 
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Percentage of binary configurations where : J(t) ■ J(0) < 1 — ej 
(2O + 1O)M0 (15 + 15)Mq (20 + 5)Mq (10 + 10)Mq {7 + 5)Mq (10 + 1.4)Mq [NS-BH] 


ej = 0.05 

17.5% 

6.0% 

33.5% 

7.0% 

3.5% 

0.0% 

ej =0.10 

2.5% 

0.0% 

11.0% 

0.0% 

0.0% 

0.0% 


TABLE III: Deviation of the total angular momentum J from its initial direction. This table shows the percentage of the 
binary configurations where J{t) ■ J(0) goes below 1 — ej, for the ej given in the first column. 


measure of the goodness of the fit, given by 

50Hz</<250Hz ' ^ ^ ^ ' 

is shown in Tab. m The maximum deviations are all smaller than ~ 0.1 rad, except for the lighter (7+5)M0 BBH and 
(10 + 1.4)M0 NS-BH systems (where however the average deviations are still ~ 0.1 rad). This suggests that templates 
with phasing expressions similar to 15211 (such as those proposed in BCVl) could already approximate rather well the 
full target model studied in this paper. 


E. Simple and transitional precession of total angular momentum 

For most of the binary configurations investigated, we find, in analogy with the ACST analysis, that the direction 
of total angular momentum does not change much during evolution. In other words, transitional precession does not 
occur. Table HTTI shows the fraction of configurations that yield 

mm J(t) • Jo < 1 — ej , (54) 

when ej = 0.05 and 0.10. Let us now try to understand the numbers of Tab. nm in more detail. 

We first focus on the columns two to six, which deal with binaries of maximally spinning BHs. For BBHs with 
single masses m = 5-2Q Mq, the total spin is not usually large enough to satisfy the transitional-precession condition 
(El), as we can prove easily by using all the evolution equations at the leading PN order: during the evolution, the 
magnitude of the orbital angular momentum decreases with the GW frequency /, as in 

|L|«|LN| = r;(7rM/)-i/3M2, (55) 


while the total spin is bounded by 


|S| < |Si| + IS 2 I =ml + ml = il- 2r])M^ . (56) 

In order for transitional precession to occur, we need at the very least |Ln| = |S| [see Eq. EJ]: which requires 

7?(^M/)-i/3 < (1 _ 2r ,), (57) 


or 


J ^ /ti 


7rM(l — 2r])^ 


(58) 


For transitional precession to occur before we reach the Schwarzschild ISCO frequency /schw = I/v^ttM, we then 
need 


/schw y 1 - 2?7 y ~ 


rj > 0.22. 


(59) 


Although the ending frequencies obtained within our target model are usually higher than /schw) the very configu¬ 
rations that can have transitional precession (those with nearly antialigned total spin and orbital angular momenta) 
have always lower ending frequencies, making 0.22 too large an estimate for the critical value of 77 . 
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FIG. 8 : Transitional precession. Evolution of the direction of total angular momentum (left panel) and of Newtonian orbital 
angular momentum (right panel) in the transitionally precessing (2O + 5)M0 BBH with initial angles 9si = 175.4°, 6 s 2 = 105.4°, 
and — (f>s 2 = 92.0° (at /gw = 30Hz). 


As a consequence, among all the configurations we have considered, only (20 + 5)Mq and (20 + 1O)M0 BBHs can 
then have observable transitional-precession phases. These latter binaries are characterized by significantly larger 
changes in J [see Tab. IIIII . However, (20 -I- 1O)M0 BBHs still require / > /™aS = 138 Hz, which is very close to the 
relevant ending frequency; so the change in J is smaller, and we never observed episodes of transitional precession 
in the 200 initial configurations analyzed. On the contrary, we observed a few for (20 -I- 5)M0 BBHs; one example 
follows from the initial configuration given by = 175.4°, 9s2 = 105.4°, and (j)Si — 4'S2 = 92.0° (at /gw = 30Hz). 
In this configuration the initial spin of the more massive body is almost exactly antialigned with the orbital angular 
momentum. The trajectories of J and Ln during this evolution are shown, respectively, in the left and right panels 
of Fig.lHl 

By contrast, none of the NS-BH configurations examined exhibits transitional precessions. This is because the BH 
is taken as maximally spinning, so S is always much larger than L in the frequency band under consideration. 


F. Apostolatos’ power law for orbital precession 


As discussed in the previous section, the vast majority of binary configurations undergoes simple precession, where 
J remains constant, while Ln and Si ,2 precess around it. For ACST configurations (toi « m 2 and negligible SS 
interactions, or S 2 ~ 0), both Ln and S precess around J with the precession frequency 0, Eq. (42)] 


Hp = 


da-p 

dt 





(60) 


ACST identified two regimes where the evolution of ap can be approximated very well by a power law in uj (or /). 
For Ln ^ S, the total angular-momentum J « Ln ~ using to ~ it is straightforward to derive from 

Eq. 16011 that Qfp is approximated well by a linear function of / ^, 

«^^_i)(/)« j+62, ( 61 ) 

where Bi and B 2 are constant coefficients. Since Ln/S ~ the condition Ln ^ S corresponds to comparable- 

mass binaries {rj ~ 1/4) or to large separations. For Ln «C S', we have J « S; in this case we derive from Eq. iCT) 
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FIG. 9: Simple precession. The upper graphs show the evolution of the direction of total angular momentum J (left), and of 
Newtonian orbital angular momentum Ln (right), in the case of the simply precessing (20 + 5)Mq BBH with initial angles 
Osi = 44.6°, 0 s2 = 101.0°, and <j)Si — 4>S2 = —39.7° (at 30Hz). The lower graphs show the projection of Ln onto the plane 
perpendicular to the initial J (left), and the angle between Ln and J, plotted as a function of inverse GW frequency (right). 
The BBH was rotated in space so that the initial direction of J would be parallel to the z axis. 


that Qfp is approximated well by a linear function of / 


«P^2/3)(/) 


p/3 


-B' 


(62) 


where B[ and B2 are constant coefficients. The condition S corresponds to toi <C m 2 or to small separations 

(late inspiral). 

It turns out that Eqs. and (P|l apply also to a large fraction of the BBHs and NS-BH binaries studied in 
this paper. This can be tested semiquantitatively by the following procedure. For each configuration, we take the 
precession angle ap{f) and we fit it with a function - 2 / 3 )(/) form (IHTll or for frequencies in the 

range 50-250 Hz. We then evaluate the maximum difference 


^Q^max(-l,-2/3) = max 

^ ' 50Hz</<250Hz 


®p(/) Q;p(-i,-2/3) (/) 


(63) 
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90% percentiles of error in precession angle, Aomax 
(15 + 15)Mq (20 + 10)Mq (20 + 5)Mq (10 + 10)Mq (7 + 5)Mq (10 + 1.4)Me [NS-BH] 


A ^ 90% 
^'^max(-l) 

0.30 

0.24 

0.23 

0.34 

0.64 

0.61 

a^90% 

^^maxf-2/3) 

0.52 

0.48 

0.50 

0.68 

1.14 

0.72 


TABLE IV: Approximation of binary precession histories using best-fit parameters Bi and B 2 in Eqs. lICTl and 1621 . This table 
shows the 90% percentiles of Aaniax(-i) [Eq. lOdll l and Aaiiiax(- 2 / 3 ) in the BBH and NS-BH populations studied throughout 
this section. 


In Tab. IIVI we show the values of (that is, the 90% percentile of Aamax(-i)) and ^°r 

(15 -I- 15)Mq, (20 -I- 1O)M0, (10 -I- 10)Mq, and (7 -b 5)Mq BBHs, and for (10 -b 1.4)M0 NS-BH binaries. The numbers 
show that Eqs. inu and EH) yield (roughly) comparable approximation. This result is confirmed also by the more 
detailed analyses discussed later in this paper. 

Figure El plots the 2PN evolutions of J (upper left panel) and Ln (upper right panel) for a (20 -b 5)Mq BBH with 
initial conditions 9si = 44.6°, 632 = 101.0°, and (j)Si — 4 >S 2 = —39.7° (at 30Hz). The figure plots also the projection 
of Ln onto the plane perpendicular to the initial J (lower left panel), and the precession angle ap between Ln and J, 
plotted as a function of inverse GW frequency f~^ (lower right panel), and showing a very nearly linear dependence. 

Building on the results obtained by ACST, Apostolatos [ij conjectured (quite reasonably) that orbital precession 
will modulate the gravitational waveforms with functional dependencies given by Eqs. ED and EUl- On the basis 
of this conjecture and of the observation that, in matched-filtering techniques, matching the phase of signals is more 
important than matching their amplitudes, Apostolatos proposed a family of detection templates El obtained by 
modifying the phasing of nonspinning PN templates as in 

Apostolatos ansatz: V^spinning ^ V^non spinning “b C COs((5 -b Bf ^ ) , (64) 

while keeping a Newtonian amplitude Recently, Grandclement, Kalogera and Vecchio da applied Apostolatos’ 

suggestion to an approximated analytical model of NS-BH binaries and low-mass BBHs: whereas the addition of phase 
modulations according to Eq. did increase the effectualness Q of nonspinning PN templates, the resulting DTE 
family was still not good enough to recommend its application when trying to capture the real modulated waveforms. 
Moreover, this DTE requires three additional intrinsic parameters (C, 5, and B) on top of the two BH (or NS) masses. 
The resulting GW searches would then be plagued by an extremely high computational cost. 

In the rest of this paper, we shall propose a better template family, inspired by old and new insight on the 
precessional effects that appear in the gravitational waveforms. As we shall see, Apostolatos’ ansatz can be improved 
to build DTFs that have both high effectualness Q and low computational requirements. 


IV. DEFINITION OF MODULATED DTFS FOR PRECESSING BINARIES 

We are now going to bring together all the observations reported in Sec. cni to build DTFs that perform well in 
capturing the detector response to the GWs emitted by precessing binaries of NSs and spinning BHs (at least as long 
as the actual physical signals are modeled faithfully enough by the adiabatic target model described in Sec. CH). 

In Sec. irai we develop a new (as far as we know) convention for the generation and propagation of GW from 
spinning binaries; this convention has the desirable property of factorizing the waveform into a carrier signal whose 
phase is essentially the accumulated orbital phase of the binary, and a modulated amplitude term which is sensitive to 
the precession of the orbital plane. In Sec. nvH we then use the results of Sec. nTml to build an approximation of the 
carrier signal, and the results of Secs. mrsiinrEi and nmi to build an approximation to the modulated amplitude; 
using these terms together, we define three families of detection templates. In Sec. ITVT1 we describe two standard 
families of nonspinning-binary templates; in Sec. I VII we shall compare their performance with the performance of our 
DTFs, to evaluate the performance improvements brought about by our treatment of precession. 


A. A new convention for GW generation in spinning binaries 

At least two conventions are used to express the gravitational waveforms generated by binaries of spinning compact 
objects, as computed in the qua drup olar approximation [^: the AGST convention |T^, which uses a rotating reference 
frame, and the EG convention |l^ . which uses a nonrotating reference frame. We discussed the EG convention in 
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Sen, inn and we used it throughout this paper to generate gravitational waveforms from the numerical integration of 
the equations of motion of the target model. Before going to the specific conventions, we shall first sketch a generic 
procedure to write the gravitational waveform. 

In general, the unit vector along the separation vector of the binary, n(t), and the unit vector along the corresponding 
relative velocity, A(t), can be written as 

n(f) = ei(t) cos$(t) + e 2 (t) sin<i)(t), X{t) = —ei(t) sin<i)(t) + e 2 (t) cos$(t), (65) 

where ei(t), e 2 (t), and 63 ( 1 ) = LAr(t) are orthonormal vectors, and ei. 2 (t) forms a basis for the instantaneous orbital 
plane [see Fig.^; the quantity $(t) is then the orbital phase with respect to The definition of ei^ 2 {t) and of 

$(t) is not unique: an arbitrary function of time can be added to $(t), and then compensated by a time-dependent 
rotation of ei 2 (t) around L 7 v(t), leaving n(t) and X{t) unchanged. In nonspinning binaries the orbital plane (and 
therefore Lat) does not process, so the natural choice is to keep 61^2 constant. In spinning binaries LAr(0 precesses, 
and different, but nonetheless meaningful, conventions can be given for ei^ 2 {t) and $(t). Note that $(t) is not^ in 
general, the same as the accumulated orbital phase 'l'(t) = J oj{t)dt. Given a convention for ei^ 2 (t) and $(t), the 
tensor that appears in Eq. d can be written as 

Q^J = —2 ([e+J*-^ cos2$ -I- sin2<i)) ( 66 ) 

where 


e+ = ei (g) ei - 02 G 62 , ex = ei (g 62 -b 62 (g ei. 


(67) 


With the detector lying along the direction N, one goes on to define a radiation frame, formed by orthonormal vectors 
e^, By and b^ — N. The GW response is then given by 


/iresp = ([e+]*^' cos2$ + [gx]*^' sin2$) ([T+]i, F+ + [Tx]y F’x) , 

^ ^ '-V-''---' 

factor Q: quadrupole moment factor P: detector projection 


( 68 ) 


where the tensors [T+_x]ii are given by (EHll . namely 


T. = e 




> ( 




- + — '-'X '-'X '^y '^y 

and where id|_ and are given by Eq. d, namely 


Tx = e. 






T+, X — r, 


I G.x Go, 


^ +, X Ju 5 


(69) 


(70) 


with ex,y the unit vectors along the orthogonal arms of the interferometer. Again, and are not uniquely defined, 
because they can be rotated at will around N, of course changing the values of F+ and . 

AGST refer to the direction N of GW propagation, by imposing that ef^®'^(t) oc N x Ljv(t); they also set 
e^{t) oc ±N X LN(t)- Although the ACST convention has allowed some insight into the waveforms, it is rather 
inconvenient for the purpose of data analysis, because almost all the quantities that come into Eq. ( 1 ^ [ei, 2 , T+_ X, 
and id|-_ x] depend both on the time evolution of the binary and on the direction to the detector. Using the terminology 
introduced in Sec. nm and Tab.m under the ACST convention the local and directional parameters are entangled in 
a time-dependent manner. 

EG introduce the fixed source axes {ef,ey,ef} [see Sec. IllCj . and they impose that ef(f) oc ef x LAr(t) [see 
Eq. (I17ll ]. The radiation frame does not change with time [see Eqs. 121ll - 128|l ]. As a consequence, the factors Q and P 
in Eq. become disentangled: factor Q expresses the components of the quadrupole moment, which depend only on 
the evolution of the binary inside the source frame; factor P expresses the projection of the quadrupole moment onto 
the radiation frame and onto the antisymmetric mode of the detector, which depend only on the relative orientation 
between the source frame and the detector. However, for our purposes there are still two shortcomings in the EC 
convention: 


1. The EC convention defines ei, 2 (t) and <I>(t) in terms of the fixed source frame y which is quite artificial, 
because only the relative orientation between binary and detector affects the detector response /iresp- 

2. In Sec. IIIIl3l we saw that the accumulated orbital phase '!'(<) is (almost) nonmodulated, so the modulations of 
the waveform come mainly from the precession of the orbital plane. Under the EG convention, the modulations 
appear only in factor Q of Eq. but they appear both in the phase $(t) and in the precession of the tensors 
e+^x( 0 - It would be nice to isolate the precessional effects in either element. 
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convention 

factor P 

T+,x F'+,x 

factor Q 

m 

e+,x(f) 

ACST 

function of basic, local, and 
directional parameters; 
time dependent 

function of basic, local, and 
directional parameters 

function of basic, local, and 
directional parameters 

FC 

function of directional parameters; 
time independent 

function of basic, local, and 
directional parameters 

function of basic, local, and 
directional parameters 

precessing 

function of directional parameters; 
time independent 

function of basic and local 
parameters only; coincides with 4'(t) 

function of basic and local 
parameters only 


TABLE V: Parametric dependence of the building elements of the detector response function /iresp [Eq. I 68 I 1 I under the ACST, 
FC, and precessing conventions. 


Both issues would be solved if we could find a modification of the FC convention where $ coincides with the accu¬ 
mulated orbital phase, T'. As it turns out, it is possible to do so: we need to redefine the vectors ei^ 2(0 so that they 
precess alongside L^V) 


ei{t) = Cleit) X ei{t), i = l,2, 


(71) 


with 


Oe(t) = — [flLit) ■ Lw(i)] LAr(t) , (72) 

where is obtained by collecting the terms that (cross-product) multiply Ljv in Eq. 0 . In App. m we prove that 
this convention yields $ = w = 'F, as desired. Qualitatively, one can reason as follows. The angular velocity of the 
binary lies along Ljv(t), and has magnitude 'F = w. The reason why $ and ill differ is that the orbital basis 61 ^ 2 , 
used to define $, must rotate to keep up with the precession of the orbital plane. However, the difference vanishes if 
we constrain the angular velocity of 61^2 to be orthogonal to L^v; Eq. 17211 provides just the right constraint. In the 
following, we shall refer to our new convention as the precessing convention. 

In Tab.0we summarize the parameter dependence of the terms that make up the detector response function [Eq. 
(P|l ]. under the three conventions. It is important to remark that in the precessing convention the polarization tensors 
e+^ xit), as geometric objects, do not depend on the source frame, but only on the basic and local parameters. In 
practice, however, we need to introduce an arbitrary choice of the source frame to relate the orientation of the binary 
to the direction and orientation of the detector (that is, to write explicitly the products [e+, x]ii[T+, x]ij)- We can 
avoid this arbitrariness by setting the source frame according to the initial configuration of the binary at a fiducial 
orbital frequency; for example, we can impose (without loss of generality) 

ef cxSi(0)-[Si(0)-Ln(0)]Ln( 0), e^=LN(0)xef, ef = Ln(0) , (73) 

and 

ei(0)=ef, 02(0) = e^, 63 ( 0 ) =ef. (74) 

[If Si(0) and Ln(0) are parallel, ef can be chosen to lie in any direction within the plane orthogonal to Ln(0).] Then 
the initial conditions, as expressed by their components with respect to the source frame, are determined only by the 
local parameters. 


Ln(0) = (0,0,1), (75) 

Si(0) = (sindsi, 0, cos 6 >sJ , (76) 

82 ( 0 ) = (sincos(<()s 2 -(()si), sinsin((()s 2 -<(>Si), cos 6 »sJ, (77) 

along with an initial orbital phase T'q given by 

n(0) = ei(0) cos'I'o -I- 02 ( 0 ) sinT'o . (78) 


With this choice, all the directional parameters are isolated in factor P of Eq. EHl), while the basic and local parameters 
(which affect the dynamics of the binary) are isolated in factor Q. We will call upon this property of the precessing 
convention in Sec. IVI PI where we propose a new family of templates for NS-BH binaries built by writing a set of 
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orthonormal component templates that contain all the dynamical information expressed by factor Q, and then using 
their linear combinations to reproduce the projection operation expressed by factor P. 

Going back to the main thrust of this section, we obtain the detector response /iresp by setting the direction to 
the detector N (specified by the angles 0 and ip with respect to the source frame), and by introducing the radiation 
frame, oriented along the axes 

ef = -ef sin(/? + e^ cos(/j, (79) 

By = — ef COS0 cos(/3 — COS 0 sin(/j + ef sin© , (80) 

ef = +ef sin 0 cos (/? + sin 0 sin (/? + ef cos 0 = N ; (81) 

we then get 

hresp = ([e+]*'' cos2^ + [ex]*^' sin2^') ([T+],, F+ + [Tx]y F^) . (82) 

Applying the stationary-phase approximation (SPA) at the leading order, we can write the Fourier transform of /iresp 
as 

kespif) =-hcif) +i[e^itf)]F) ([T+],fc F++ [Tx],fc ) for / > 0 , (83) 

where hcif) is the SPA Fourier transform of the carrier signal, 

he = — cos 2'ii, (84) 

U r 

and where t/ is the time at which the carrier signal has instantaneous frequency /. 


B. Definition of a new DTF for processing binaries 


By adopting the precessing convention, we isolate all the modulational effects due to precession in the evolving 
polarization tensors (these effects will show up both in the amplitude and in the phase of /iresp)- The 

discussion of Sec. IIII PI shows that, to a very good approximation, the carrier signal is not modulated, so we expect 
that hcif) should be approximated well by the nonspinning PN templates studied in BCVl, or variations thereof. 
As for the time dependence of the tensors [e+^x]*'^ , the discussion of Secs. IIII El and HITU suggests that we adopt the 
Apostolatos’ ansatz |^, and write expressions in the generic forms 

[e+,x]*^[T+,x]jfc ocC+,xCos(;B/“^/^ + (5+,x) or oc C+,x cos + (5+,x) • (85) 

Indeed, our extended numerical investigations provide evidence that expressions of the form (IS31) should work quite 
well for the binaries under consideration. 

All these elements suggest that we introduce a family of detection templates of the general (Fourier-domain) form 


h{ljj-NM,Ak,to, Oik', /) 


-h iak+n)Ak(f) 


g27ri/togii/’NM(/) 


(for / > 0) 


( 86 ) 


[and /i(/) = h*i—f) for / < 0], where the Akif) are real amplitude functions, the ak are their (real) coefficients, and 
to is the time of arrival of the GW signals. The function '0 nm represents the phase of the unmodulated carrier signal; 
we write it as a series in the powers of 

V'nm(/) = (V'O + ^1/2/'/' + + V'3/2/ + ...)■ ( 87 ) 

As discussed in BCVl, this phasing works well for relatively high-mass, nonspinning BBHs, and for NS-BH binaries; 
in addition, as anticipated in Sec. ITTT PI the PN coefficients are able to capture the nonmodulational effects of 
spin-orbit and spin-spin couplings on the orbital phase. In this paper we examine three specific families of detection 
templates of this form, listed in Tab. ED The subscripts “2”, “4”, and “6” in our abbreviations for the template 
families denote the number of ak coefficients that appear in Eq. 

The families (V’o'*/' 3 / 2)2 and ('0o'03/2Q!)4 were already studied in Ref. [i| for the case of nonspinning binaries. Both 
families contain the leading f~'^A Newtonian dependence of the amplitude; however, (V’o'*/' 3 / 2 <a )4 contains a correction 











21 


Template family 

V'nm)/) 

Al(f) 

A2if) 

Asif) 

{i>0lp3/2)2 

(■/oV’3/2a)4 

(/'oV’3/2S)6 

(/'0Z^3/2S')6 

zA0/-®/^+z/3/2/-^/^ 

/-"/®e(/cut-/) 
/-"/®e(/cut-/) 
eifeut - f) 
/-"/®e(/cut-/) 

eifeut - f) 

/-"/® COs(I3/-2/3)0(/cut -/) 
cos(I3/-i)0(/cut-/) 

/-^/«sin(S/-2/3)0(/,,t-/) 

sin(B/-i)e(/cut-/) 


TABLE VI: Specification of the DTFs examined in this paper. 


to the Newtonian amplitude (introduced in BCVl, where it was parametrized by a) which can account for the variation 
of the rate of inspiral in the late stages of orbital evolution. The first family is given by 

(^oV' 3 / 2)2 : h{...-J) = (ai +za 2 )/"’^/® 6 '(/cut - /) exp f + z/> 3 / 2 /“^^^]; (88) 

here ai +ia 2 can also be written as ^exp , where (f>2^ is the initial GW phase, and A is an overall normalization 
factor for the template. So the two at coefficients encode the initial global phase of the waveform, plus a normalization 
factor. The second family is given by 

(' 0 o' 03 / 2 a )4 : h{...;f) = [(oi +102)/"’^^® + (03 + *<24)/"^^^] 6 '(/cut - /) exp + V'3/2/”^^^]; ( 89 ) 


another way to rewrite the coefficients 0 : 1-4 more physically is ylexp[z())(f'^]/“^/®(l + oexp[z^“]/^/^), where a is the 
additional amplitude parameter and ((>“ is the relative phase of the amplitude correction (as in BCVl, in this paper 
we always set ())“ = 0). So the four coefficients Ofc encode the global phase, the strength of the correction to the 
Newtonian amplitude, and the relative phase of this correction with respect to the Newtonian amplitude, plus an 
overall normalization factor. 

The third family, (z/’o'03/2^)6j contains the leading Newtonian amplitude, modified by two modulation terms [a 
generalization of the Apostolatos’ ansatz ®] that account for the precession of the orbital angular momentum due 
to spin effects. It is given by 

(V'o' 03 / 2^)6 : = /“'^/®[(q:i+ 102) + (as+ zq;4)cos(;B/“^/^) + (05+ za6)sin(;B/“^/^)] 

X 0{fcut - /) expi[V:o/-®/" + z/^3/2/-'/"]; (90) 

another way to rewrite the six coefficients ai-g in close analogy to Apostolatos’ ansatz is 




1 + C cos{Pf + i5i + iS 2 ) 

1 + Ceos cos(/ 3/-"/3) + sin(/3/-2/3) 


(91) 


(where all the coefficients are still real). So the six coefficients Uk encode the global phase, the strength of the 
amplitude modulation, its relative phase with respect to the Newtonian amplitude, and the internal (complex) phase 
of the modulation. It is clear that our family implements a generalization of Apostolatos’ ansatz, because we allow 
a complex phase offset between the Newtonian and the sinusoidal amplitude terms, and also between the cosine and 
sine modulational terms. We consider also a variant ('0o'03/2‘BO6 of this family where the frequency dependence 

in the sinusoidal amplitude functions is replaced by f~^. 

For all three families, the templates are terminated at a cut frequency /cut) above which the amplitude drops to 
zero; this /cut is in effect one of the (intrinsic) search parameters. For all three families, the frequency dependence of 
the phase includes the leading Newtonian term, and a term that corresponds to the 1.5PN correction 

in the phase evolution of nonspinning binaries (as obtained, in the SPA, by integrating the energy-balance equation 
through an adiabatic sequence of circular orbits, using PN expanded energy and flux). In BCVl we found that 
including either the IPN or 1.5PN term is in general sufficient to model the phase evolution of nonspinning binaries 
of high mass. 


C. Definition of the standard SPA template families 

In this section we define two families of standard nonspinning-binary templates, obtained by solving the Taylor- 
expanded energy-balance equation for an adiabatic sequence of quasicircular orbits, and using the stationary-phase 
approximation (SPA) to express the result as a function of the GW frequency / (see BCVl). In Sec. IVIl we compare 
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the matching performance of these templates to the performance of our new DTFs, to show that the various tricks 
used to build the new families do indeed improve their effectualness . The standard SPA families are built from the 
analytic expressions of Refs. miii The frequency-domain phasing (under the assumption of nonevolving orbital 
angular momentum and spins) is given by El 


o [20 /Tdt 11 \ 

V’Spa(/) = 2TTftc-(l)c + — (ttM 1 -f — \^— + —r]j {ttM /)2/3 _ 4(47r - SO) (ttM /) 


/3058673 5429 617 , \ .^4/3 

+ 10 ( + 777^ - SS (ttM /)4/3 


V 1016064 1008 


144 


(92) 


where A4 = is the chirp mass, and where SO and SS are the spin-orbit and spin-spin terms, given explicitly 

by 


SO 

SS 


1 

M2 


/113 25 m2\g /113 25 mi\ 



1 

48toi 7712 M2 


-247 Si • S 2 + 721(Si • In) (S 2 • Ln) 


(93) 

(94) 


We neglect all PN corrections to the amplitude, by adopting its Newtonian functional form, we also neglect 

all precessional effects, by setting SO = SS = 0. Templates of this form are routinely used in searches for GW signals 
from nonspinning binaries. In that case, the templates are generally ended at the GW frequency corresponding to 
the Schwarzschild ISCO /schw — 0.022/M. We denote such templates as SPAs. We introduce also a variant of this 
family, SPAc, characterized by the additional frequency-cut parameter /cut, used also in our DTFs. Altogether, we 
get 

SPAs : h{M,r],to,'ij}o,aN; f) = Q!Ar/“’^/®6»(/schw -exp7[V’sPA + V'o]; (95) 


SPAc : h{M,r], fcnt,to,'ipo,aN; f) = expi[^sPA + V'o]- (96) 


V. GW DATA ANALYSIS WITH THE DTE 

In searching for GW signals using matched-filtering techniques, we construct a discrete bank of templates that 
represent all the possible signals that we expect to receive from a given class of sources. We then proceed to compare 
each stretch of detector output with each of the templates, computing their overlap (essentially, a weighted correlation). 
A high value of the overlap statistic for a given stretch of detector output and for a particular template implies that 
there is a high probability that during that time the detector actually received a GW signal similar to the template. 
This technique is intrinsically probabilistic because, for any template, detector noise alone can (rarely) yield high 
values of the statistic. In general, the higher the value of the statistic, the harder it is to obtain it from noise alone. 
So it is important to set the detection threshold (above which we confidently claim a detection) by considering the 
resulting probability of the false alarms caused by noise. 

To verify whether the DTFs developed in Sec. m can be used to search reliably and effectually for the GWs from 
spinning binaries, we need to evaluate the fitting factor FF of the DTFs in matching the target signals for a variety 
of binary and detector parameters. The FF is defined as the ratio between the overlap of the target signal with the 
best possible template in the family and the overlap of the target signal with itself |^. So in Sec. EH we discuss 
the maximization of the overlap over template parameters for a given target signal. The other important element 
to evaluate the reliability and effectualness Q of the DTFs are the detection thresholds that the DTFs yield for a 
given false-alarm probability. In Sec. EH we discuss these thresholds under the simplifying hypothesis of Gaussian 
detector noise. The material presented in this section builds on the treatment of matched-filtering data analysis for 
GW sources given in Sec. II of BGVl (which is built on Refs. mill) , and it uses the same notations. 


A. Maximization of the overlap over template parameters 

Among all the template parameters that appear in Eq. we are going to treat the ifi, /cut and B as intrinsic 
parameters; and the ak and to as extrinsic parameters: that is, when we look within one of our DTFs for the template 
that best matches a given target signal, we will need to consider explicitly many different values of the t/j, of /cut, 
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and of B', however, for any choice of these parameters, the best ak and to are determined automatically by simple 
algebraic expressions (see Sec. II B of BCVI). For the next few paragraphs, where we discuss the optimization of the 
coefficients afc, we shall not indicate the dependence of the templates on the intrinsic parameters. 

For a given signal s, we seek the maximum of the overlap, 

max(s,/i(to,afc)), (97) 


under the normalization condition 


{h{to,ak),h{to,ak)) = 1 


(98) 


[this condition is necessary to set a scale for the statistic distribution of the overlap between a given template and 
pure noise]. Here the inner product {g, h) of two real signals with Fourier transforms g, h is defined by 


{g, h) = 2 



Su{\f\) 


d/ = 4 Re 



Sn{f) 


(99) 


(see BCVI). We proceed constructively: first, we build a new set of amplitude functions Ak{f) that are linear 
combinations of the Akif), and that satisfy the orthonormality condition (^i(/), Aj{f)) = Sij for f, j = 1, 2,... n; we 
then define an orthonormal set of single-^fc templates, 

hk{to-,f) = , hk+n{to; /) = f (for / > 0) (100) 

[and hk{f) = for / < 0], which satisfy {hi{to), hj{to)) = Sij (with z, j = 1, 2, ..., 2n) for any to. The maximized 

overlap [Eq. ll?7fll ] can now be rewritten as 

2n 

max (s,/i(to, afc)) = maxmaxY^ ak{s,hk{to)), (101) 

to,ock to ak ^' 

k—1 

while the condition 10811 is now simply “ 1- The inner maximum of Eq. IIOIII (over the ak) is achieved when 


and the maximum overlap itself is 


d/c 


(s, hk (ffi)) 


( 102 ) 


2 n 


max (s, h{to, ak)) 


= maxmax> ak{s,hk{to)) 

to &k 

k—1 


\ 


2n 

inax^(s,/ij(to))" 


(103) 


This happens essentially because the sum in Eq. itTUTll can be seen as a scalar product in a 2n-dimensional Euclidean 
space, which is maximized when the unit 2n-vector ak lies along the direction of the 2n-vector {s,hk{t)). The 
quantities (s, hj{to)) for j = 1, 2, 3,..., n are given by the two related Fourier integrals 


{s,hj) = 2Re 

(sj = 2 Im 


f* + 00 


Mf) 




{/) 


Sh{f) 




AAf)A 


A*if) 


,2mfU 


Shif) 


df. 


(104) 

(105) 


We now go back to discussing the full set of template parameters. The relevant measure of the effectualness 0 of 
a template family at matching a physical signal s is the fitting factor FF, 


FF = max 

^0 jOik if cut ,'lpi 


{s,h{to,ak)) 


(106) 


(see, for instance. Sec. II of BCVI) which is maximized over the ak, but also over the time of arrival to (also an 
extrinsic parameter), and over all the intrinsic parameters, ifi, font, and B. The fitting factor is a function of the 
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physical parameters of the physical signal s, and of course of the template family used to match it. We define also 
the signal amplitude SA for a given signal, 

SA = (107) 

SA gives the optimal overlap obtained for a template that is exactly equal to the signal (except for its normalization), 
and it is inversely proportional to the luminosity distance to the source; where we do not indicate otherwise, we always 
assume the fiducial distance do = 100 Mpc. 

The maximization of the overlap over to is easy to obtain, because the integrals (irnni and I105II can be evaluated at 
the same time for all the to using Fast Fourier Transform techniques |32| . On the other hand, the maximization over 
/cut and over the other intrinsic parameters is obtained by an explicit search over a multidimensional parameter range, 
where we look for the maximum of the partially maximized (over extrinsic parameters) overlap, given by Eq. JHIJI. 
For all the actual searches discussed in this paper we employ with good results the simplicial algorithm amoeba (3,111 . 


B. False-alarm statistics of the DTFs 


In the practice of GW data analysis, template families are used to build discrete template banks parametrized by a 
discrete set of ntuples of the intrinsic parameters. Then each of the templates is correlated with the detector output, 
to see if the detection statistic [in our case, the partially maximized correlation E3] is greater than the detection 
threshold. It is important to notice that the statistic is already maximized with respect to the extrinsic parameters, 
while the intrinsic parameters serve as labels for each of the templates. Therefore, we are effectively setting up a 
separate detection test for each of the templates in the bank. 

In this section we are going to evaluate the false-alarm probability for one such test, defined as the probability that 
detector noise alone will yield an overlap greater than the detection threshold. The total false-alarm probability is 
then obtained by multiplying the false-alarm probability for a single template by the number Afshapes of independent 
signal shapes (generally of the same order of magnitude as the number of templates in the bank), and by the number 
A/bmes of possible times of arrival to , distanced in such a way that the displaced templates are essentially orthogonal 
[ 3 ^. At the end of this exercise, we are going to set the detection threshold so that the total false-alarm probability 
is acceptably low. 

Under the assumption of Gaussian noise, the inner product {n,hj) of noise n alone with a normalized template 
component hi is (by construction) a Gaussian random variable with zero mean and unit variance (see, for instance. 
Sec. II of BGVI). Because (for the same to and for the same intrinsic parameters) all the hj are orthogonal, the 
inner products {n,hj) (for j = l,...,2n) are all independent normal variables. It follows that the statistic X = 
maxto.ofc {n, h(to, Ofe)) [see Eq. 110311 ]. given by the square root of the sum of their squares, follows the x distribution 
with 2n degrees of freedom, characterized by the probability density function and cumulative distribution function 


PDFx(2n)(-^ = x) 


^2n-lg-a;^/2 

2"-ir(n) 


GDF^(2„)(X<x) 


r(n, 0, x'^/2) 

TW 


(108) 


where we have used the generalized incomplete gamma function F(n, zo, zi) = t'^~^e~*dt. For n = 1 we obtain the 

Rayleigh distribution, typical of the maximization of the amplitude of signals with two quadratures. 

In Tab. IVIll we show the thresholds needed to obtain a total false alarm probability of 10“^, with A/times = 3 10^° 
(typical of about three years of observation with LIGO), and with the AAshapes given in the first column. We observe 
that each time we increase A4hapes by one order of magnitude, the threshold increases by about 2% (this happens 
uniformly for all n’s). On the other hand, each step in n increases the threshold by about 4%. Thus, when we design 
DTFs we should keep in mind that the best possible overlap increases with the number of templates employed, and 
with the complexity of the templates (clearly, the complexity of our DTFs increases with the number of amplitude 
functions); but the detection threshold increases as well, reducing the number of signals that pass the detection test. 
So in principle we are justified in using more numerous and more complex templates only if the gain in the overlap is 
larger than the increase in the detection threshold. 

The prospects shown in Tab. rvTTI for the models with n = 2 and n = 3 improve somewhat if we constrain the 
values that the ak can attain when they are (algebraically) maximized. We can do this, for instance, if we judge that 
certain combinations of the Ofc correspond to unphysical waveforms, but then we must be consistent and exclude any 
detections that cross the threshold within the excluded parameter region. In any case, we should remember that our 
study of false-alarm statistics is based on the idealization of Gaussian noise, which will not be realized in practice: 
real-world data-analysis schemes relie on matched-filtering techniques complemented by vetoing schemes |3^ . which 
remove detection candidates using nonlinear tests on the signal. Therefore, any DTF should be evaluated in that 
context before it is excluded for producing excessive detection thresholds within the Gaussian analysis. 
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*^^hapes 

Threshold for false-alarm probability = 10 ® 


i'fo ^'3/2)2 

(V’O ^>3/2 0)4 

(V'O ’^3/2 S)6 


n = 1 

n = 2 

n = 3 

10^ 

8.44 

8.87 

9.22 

10® 

8.71 

9.13 

9.48 


8.97 

9.39 

9.73 

10® 

9.22 

9.63 

9.97 

10® 

9.47 

9.87 

10.21 


TABLE VII: Detection thresholds for a false-alarm probability = 10“® for a x-distributed detection statistic with 2n degrees 
of freedom, for Atimes = 3 10^°, and for the Afshapes given in the first column. The values given for (V’o, V’3/2 01)4 do not take 
into account the 0“ = 0 constraint. 


VI. EVALUATION OF DTE PERFORMANCE 

We wish to investigate the effectualness 0 of our DTFs in matching the GW signals generated by precessing 
binaries of spinning compact objects, at least as approximated by the target model described in Sec. m To do so, 
we shall evaluate the fitting factor FF [Eq. HI Ofijl ] of the DTFs over a population of binaries with a variety of basic, 
local and directional parameters, and compare the results with the FF obtained for the standard SPA families [Sec. 
iTvn . In Sec. Em we study the effect of the directional parameters on FF (and SA), with the aim of reducing 
the dimensionality of the test populations. In Sec. Em we describe the Monte Carlo scheme used to generate the 
populations, and we identify two performance indices for the template families (namely, the simple and SA-weighted 
averages of FF). In Sec. IVK^I we give our results for these indices, focusing first on the BBHs considered in this paper. 
Finally, in Sec. Em we give our results for NS-BH binaries, and we briefly describe a new, very promising family of 
templates for these systems, suggested by the insights accreted during the development of this paper. 


A. Effect of directional parameters on FF and SA 

As we have seen in Secs. ESI and nroi the detector response hresp is a function not only of the basic and local 
parameters of the binary (which describe respectively the masses and spin magnitudes, and the initial relative direc¬ 
tions of the spins and the orbital angular momentum, and therefore change the dynamical evolution of the binary), 
but also of the directional parameters (which describe the relative direction and orientation of binary and detector, 
and alter the presentation of the precessing orbital plane of the binary with respect to the direction and orientation 
of the detector). Thus, all the parameters will affect both the amplitude SA = {h^esp, hresp)^^^ of the signals received 
at the detector and the ability of our DTFs to match them, as codified in the fitting factor FF; it is therefore clear 
that, in evaluating the effectualness of our DTFs at matching the target signals, we will need to compute FF not only 
for a range of binary masses and spins, but also for a suitable sampling of the local and directional parameters. 

In the case of nonspinning binaries (see BCVl), there are no local parameters as we defined them in this paper; 
the directional parameters do change the GW signal, but only by multiplying its amplitude by a constant factor, and 
by adding a constant offset to its phase (as opposed to modulating amplitude and phase as in the case of spinning 
binaries). In BGVl (following a common practice in the GW data-analysis literature), we included the variation 
of the amplitude in the definition of the target signals, by averaging the amplitude factor over uniform solid-angle 
distributions of the directional parameters [see Eq. (29) of Sec. II D]. As for the initial phase of the signal, we defined 
the FF on the basis of minmax overlaps |3|, which are maximized over the initial template phase (and over all the 
other extrinsic and intrinsic template parameters) but minimized over the initial signal phase; this minimization is 
obtained algebraically, just as for extrinsic template parameters. In fact, it turns out that minimizing or maximizing 
the overlap over the initial signal phase changes the resulting FF by a very small quantity. 

In the case of the spinning binaries examined in this paper, this picture changes radically, because minimizing 
the overlap over the directional parameters yields very low FFs that are not representative of the typical results 
that we would get in actual observations. So we take a different approach: we study the distribution of FF for a 
population of binaries characterized by different basic, local, and directional parameters. In particular, we select 
several astrophysically relevant combinations of basic parameters, and we sample randomly (but as uniformly as 
possible) the space spanned by the local and directional parameters. In practice, we can exploit certain symmetries of 
this space (that is, the fact that different combinations of the local and directional parameters yield the same signal) 
to reduce its effective dimensionality. Let us see how. 
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Under the FC convention, the complete specification of a target signal requires (at least formally) 15 parameters: 

according to our classification ISec. lllL)ll . four of these are the basic parameters (M, 77 , ^i, and S 2 ); three are the 

local binary angles {0Sn (pSi — <(> 52 ); three are the directional binary angles (0 Lni 4>Si + <(’ 52 ); 

five are the directional GW and detector angles (0, (p, 6, (/>, and V')- Of the latter, 6, (j), and tp come into the waveform 
only through the antenna patterns and [see Eqs. (123 and (EOll]. It is redundant to specify both the directional 
binary angles (which determine the orientation of the binary as a whole in space) and the directional GW angles 
(which determine the direction N of GW propagation to the detector), because if we apply the same rotation to N 
and to the binary vectors Lat, Si, and S 2 , we do not change the response of the detector /iresp- So we can use this 
freedom to set 0 = 7 r /2 and (/? = 0. Once this is done, we still have the freedom to rotate the detector-binary system 
around the axis N. Such a rotation (by an angle v) will transform the F^ and F^ antenna patterns according to 

F_|_ ^ cos 277 — Fx sin 2 z 7 , (109) 

Fx ^ F_|_ sin277 + Fx cos277. (110) 

Looking at Eqs. i|23 and we see that, for any original 9, (p, and t/;, we can always find an angle v for which 
F+ = 0. The corresponding new Fx becomes 

Fx = (1 + cos^ 9y cos^ 2(j) + dcos^ 0sin^ 2(j ); (HI) 

once again, the detector response does not change. Eor future use, let us define as p[Fx] (with p[Fx]dFx = 1) the 
probability density for |Fx | induced by uniform solid-angle distributions for 9 and (p [notice that ip does not appear 
in 111111 ) 1 . 

Now, for a given DTE and for given basic parameters, consider the distribution of FF and SA obtained for an 
11 -parameter population of target signals specified by uniform solid-angle distributions of 0 Ln,Si.S 2 j 0 Ln.Si,S 2 j 'Fj 
9, (p, and ^p. By the above arguments, we obtain the same distribution of FF and SA from a 6 -parameter population 
of target signals specified by uniform solid-angle distributions of 6 *Ln,Si,S 2 i 0 Ln,Si.S 2 ! by 0 = 7 r/ 2 , (/? = 0, F+ = 0, 
and by Fx distributed acoording to p[Fx]. Moreover, because Fx appears only as a normalization factor in front 
of the expression m for the signal (once F+ = 0), we can simply set Fx = 1: this operation does not change FF 
[because Fx appears homogeneously in the numerator and denominator of Eq. 11061) ]. while the distribution of SA 
for the original 11 -parameter population can be recovered from its moments on the 6 -parameter population: 

(SA™)n-pa, = / / (Fx)™SA>[Fx]dFx\ = {SA^h-p^r f (Fx)>[Fx]dFx. (112) 

/ 6 -par -^0 


B. A Monte Carlo procedure to evaluate DTF performance 

We are going to evaluate the effectualness Q of our DTFs within a Monte Carlo framework, by studying the 
distribution of FF (and FF^SA^, see below) over six sampled populations of 1000 binaries each, specified as follows. 
We study the binary systems already examined in Sec. irm BBHs with masses (20 + 10)Mo, (15-fl5)Mo, (20-f5)MQ, 
(10 -I- 1O)M0, and (7 -I- 5)Mq, and NS-BH binaries with masses (10 -I- 1.4)M0. All the BHs have maximal spin, while 
the NSs have no spin. We integrate numerically the target-model equations starting from initial configurations that 
correspond to instantaneous GW frequencies of 30 Hz when M > 20Mq, and 40 Hz otherwise. For each set of masses, 
we use the Halton sequence with bases 2, 3, 5, 7, 11, and 13 to generate 1000 quasirandom sets of the six angles 
^Ln,Si,S 2 ()'Ln,Si,S 2 ; directions of the resulting orbital angular momentum and spins are uniformly distributed 
over the solid angle. We denote each sestuple by the sequential index I, for I = 1,... ,J\f = 1000. We always set 
0 = n/2, (p = 0, F+ = 0, Fx = 1, and we take do = 100 Mpc. 

For each set of masses, and for each DTF, we compute the Monte Carlo average of the FF, 

_ 1 ^ 

FF=(FF) = ^^FF[;], (113) 

^ 1=1 


and its variance 


1 


E 


z=i 


(FF[^]-FF)^ 


= (AFF^) 


(114) 
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which can be used to estimate the sampling error of the Monte Carlo average as AFF ~ ctff/vW- 

There is another function of FF and SA that has a particular interest for our purposes. Consider each configuration 
I as a representative of a subclass of physical signals that have the same binary, GW, and detector parameters (except 
for the degenerate parameters discussed above), but that are generated uniformly throughout the universe. The rate 
of successful signal detections using a given DTF is then 

= = (115) 

where TZdo is the rate of events out to the distance do from Earth. Here we assume that TZdo is a function of the basic 
parameters of the binary, but not of 1. This equation holds because FF[Z] SA[1] is the signal-to-noise ratio (that is, the 
overlap maximized over the DTF) for the signal I at the distance do; the ratio of FF[Z] SA[1] to the DTF threshold 
gives the fraction or multiple of the distance do out to which signals of the class I will pass the detection test. Folding 
in p[Fx] we get 

F-detectll] = F-detectl^Fx = ^] ' [ (Fx)^ p[Fx]dFx = 0.293 • TT-detect [^, = !]■ (116) 

Jo 

Summing over the I, we get an estimate of the total detection rate, TT-detect = (1/Af) ^detect[^]- On the other 
hand, the optimal detection rate that we would obtain with a perfectly faithful DTF is 


"^optimal — F-do 




\ threshold [DTF] J 


{Fxfp[Fx]dFx 


(117) 


We can therefore define the effective average fitting factor FFes (which is a function of the basic parameters of the 
binary, but which is already integrated over 1) from the equation 


*77-detect — h'Fgg-TT-optiinal. 


(118) 


We then get 




(119) 


To compute the Monte Carlo results presented below we use the jackknifed version of this statistic to remove bias, 
and we estimate the error AFFefi as the jackknifed sampling variance. For each class of binaries and for a specific 
DTF, the effective fitting factor FFes represents the reduction in the detection range due to the imperfection of the 
DTF. The corresponding reduction in the detection rate is FFgfj. 

In Fig. ^1 we show two examples of the distribution of signal amplitudes for the (15 + 15)Mq BBHs and for the 
(10 + 1.4)Mq NS-BH binaries in our Monte Carlo population (as computed with the 2PN target model). The plots 
show SA as a function of the initial Jat • N, normalized at distances that yield SAs comparable to typical detection 
thresholds, and averaged over the probability distribution p[F'x]. For heavy, comparable-mass BBHs (except perhaps 
for the last stages of the inspiral), the orbital angular momentum Ltv is much larger than Si_ 2 , so the initial total 
angular momentum Jjv is almost perpendicular to the orbital plane; furthermore, as seen in Sec. IHI El the direction 
of Jat does not change much during evolution. Because in the quadrupole approximation the emission of GWs is 
stronger along the direction perpendicular to the orbital plane, values of [Jat • N| close to one give stronger signals, as 
seen in the left panel of Fig.^| For NS-BH binaries, where p is small, the BH spin Si is much larger than Lat, and 
Jat lies roughly along Si. So the upward curve of the left panel appears when Lat is roughly parallel or antiparallel 
to Si and Jat (that is, when the conserved quantity KeS oc Lat • Si has a large absolute value), while a downward 
curve appears when Lat is orthogonal to Si and Jat (that is, when HeS has a value close to zero |37jl. The mixture 
of these two tendencies creates the shape seen in the right panel of Fig. HOI 


C. Performance indices for the standard SPA templates and for the modulated DTFs 

Figure m shows the distribution of FFs, evaluated for our DTFs and for the SPA standard templates against 
the 2PN target model, within the Monte Carlo populations of BBHs and NS-BH binaries described in the previous 
section. The vertical lines show the Monte Carlo estimates of FF and FFe® (the latter is always larger), with their 
estimated errors; these numbers are given also in Tabs. IwTl andinn We wish to discuss several features of the FFs. 
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FIG. 10: GW signal amplitude SA as a function of the initial Jat • N (that is, the cosine of the angle between the direction 
of GW propagation and the initial total angular momentum at the Newtonian order), for our Monte Carlo populations of 
(15 + 15) Mq BBHs (in the left panel) and (10 +1.4)M© NS-BH binaries (in the right panel). The signal amplitude is computed 
for a LIGO-I noise curve [Eq. (28) of BCVl]; it is normalized at fiducial distances of 100 and 30 Mpc, and averaged over the 
probability distribution p[Fx]- 



Fitting factors against 2PN target model 

(7 + 5)M0 (1O + 1O)M0 (15 + 15)M0 (20 + 5)M0 (20 + 1O)M0 

FF FFeff FF FFeff FF FFeff FF FFes FF FFeff 

SPAs 

SPAc 

(V’oV’3/2)2 

('!/'0/'3/2a)4 

{'<p0tp3/2B)6 

0.9030(24) 0.9390(15) 
0.9018(23) 0.9367(18) 
0.9262(22) 0.9595(13) 
0.9288(22) 0.9617(13) 
0.9753(07) 0.9828(05) 

0.8944(21) 0.9198(12) 
0.9294(20) 0.9558(12) 
0.9423(17) 0.9657(10) 
0.9480(16) 0.9703(10) 
0.9861(03) 0.9895(02) 

0.8105(25) 0.8282(16) 
0.9313(18) 0.9548(10) 
0.9414(15) 0.9620(08) 
0.9551(14) 0.9726(08) 
0.9863(03) 0.9891(02) 

0.8576(25) 0.8844(22) 
0.8854(23) 0.9096(21) 
0.8921(22) 0.9178(23) 
0.8986(21) 0.9212(23) 
0.9746(05) 0.9794(05) 

0.8264(27) 0.8494(18) 
0.9186(20) 0.9461(12) 
0.9270(17) 0.9529(12) 
0.9421(16) 0.9625(12) 
0.9843(03) 0.9884(03) 


Fitting factors against 3.5PN target model 

(7 + 5)M0 (1O + 1O)M0 (15 + 15)M0 (20 + 5)M0 (20 + 1O)M0 

FF FFeff FF FFeff FF FFeff FF FFeff FF FFeff 

('i/>0/’3/2l3)6 

0.9708(08) 0.9802(06) 

0.9854(03) 0.9887(02) 

0.9854(03) 0.9883(03) 

0.9738(06) 0.9775(05) 

0.9844(03) 0.9882(02) 


TABLE VIII: Averages FF and FFefi of the fitting factor FF against the 2PN and 3.5PN target models, for the DTFs and for 
the standard SPA template families, as computed on our BBH Monte Carlo populations. The numbers in parentheses give the 
estimated Monte Carlo errors on the last two digits of FF and FFeff. 


1. The SPA template families (solid and long-dashed black lines) always give the worst performance. Except for 
the lighter systems, (7 + 5 )M 0 BBHs and (10 + 1 . 4 )M 0 NS-BH binaries the SPAs family (solid black line) 
is consistently less effectual than SPAc, because the target-model ending frequencies are usually different from 
the Schwarzchild-ISCO frequencies used to terminate the SPAs templates (in the majority of cases, they are 
higher). The improvement [SPAs to SPAc] in FF is ~ 3% for M ~ 2O-25M0, and ^ 10% for M = 30Mq. As 
pointed out in BCVl, it is important to add the frequency-cut parameter /cut whenever the ending frequency 
is not known very well, but it is expected to fall within the band of good interferometer sensitivity. 

2. Although the {'ipo 'ip 3 i 2)2 DTF (short-dashed green lines) is essentially a reparametrization of SPAc (both families 

have the /cut parameter), it is slightly more effectual. The reason for this is that the physical ranges of M and 
r] used to optimize FF (and in particular the constraint r] < 0.25) limit the ability of the expression ■i/sPA(/) 
to reproduce the phasing of the target. On the contrary, in the (V'o ^ 3 / 2)2 DTF the coefficients of and 

/“^/^ are not functions of M and rj, but free phenomenological parameters that can achieve the best possible 
values to match the target phasing. This added freedom does not buy a dramatic improvement for the spinning 
binaries studied in this paper, because the SPAc templates are already rather close to the adiabatic target model 
(except of course for processional modulations). On the contrary, in BCVl we saw that using unconstrained 
phenomenological parameters with extended ranges is very important to follow the nonadiabatic dynamics of 
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FIG. 11: Distribution of fitting factor FF against the 2PN target model for the DTFs and for the standard SPA template 
families, for onr BBH and NS-BH Monte Carlo populations. The vertices of the segmented curves show the number of samples 
(out of 1000) for which the FF falls within the equispaced bins [0.725, 0.75), [0.75, 0.775), ... (the bins are plotted logarithmically 
to emphasize the region of FF close to one; notice that the NS-BH figure in the bottom right corner shows a different bin 
range). The vertical lines show the averages FF and FFes with their Icr error bars (FFes is always the larger number). 


late inspiral, as predicted by some PN models for nonspinning binaries. 

3. The ('(/’o ^ 3/2 ct )4 DTF (dot-dashed blue lines) introduces the amplitude-remodeling coefficient a. In BCVl we 
found that a (together with the extension of parameter ranges) helped follow the nonadiabatic dynamics of 
some target PN models [see Tab. m In this paper, however, the only target model is obtained in the adiabatic 
limit, so the frequency-domain amplitude (except of course for the modulations due to precession) is always 
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2PN target model 
(10 + 1.4)M0 

FF FFeff 

SPAs 

SPAc 

(0003/2)2 

(0003/28)6 

0.7800(34) 0.8169(37) 
0.7747(49) 0.8129(54) 
0.7807(41) 0.8316(46) 
0.9331(15) 0.9452(14) 


3.5PN target model 

(0003/28)6 

0.9263(15) 0.9378(14) 


TABLE IX: Averages FF and FFeft of the fitting factor FF against the 2PN and 3.5PN target models, for the DTFs and for 
the standard SPA template families, as computed on the (10 + 1.4)M0 NS-BH Monte Carlo populations. The numbers in 
parentheses give the estimated Monte Carlo errors on the last two digits of FF and FFefi. 


FF against selected BCVl PN models, for the SPAc and ('!/’oV’ 3 / 2 C*) template families 
T(2,2) T(3,3.5,e = 2) P(2,2.5) P(3,3.5,0 = 2) EP(2,2.5) EP(3,3.5,0 = 2) 

SPAc ('0 oV>3/2Q)4 SPAc (V’oV)3/2Q)4 SPAc (V’oV) 3/2Q)4 SPAc (V’oV) 3/2Q)4 SPAc ('ipoi>3/2a)4 SPAc (V’oV)3/2Q)4 


(1O-MO)M0 

0.984 

0.992 

0.984 

0.988 

0.979 

0.985 

0.959 

0.990 

0.988 

0.994 

0.949 

0.994 

(20+5) Mq 

0.970 

0.992 

0.960 

0.986 

0.950 

0.978 

0.968 

0.985 

0.930 

0.993 

0.967 

0.993 

{20+10)Mq 

0.964 

0.989 

0.959 

0.986 

0.925 

0.977 

0.964 

0.986 

0.978 

0.993 

0.982 

0.993 

{15+15)Mq 

0.939 

0.989 

0.941 

0.987 

0.931 

0.980 

0.967 

0.987 

0.971 

0.991 

0.983 

0.991 


TABLE X: Fitting factors against selected PN models of nonspinning binaries (see BCVl), for the SPAc and (V’o V’ 3/2 0)4 
template families. Notice that the ('i/>o'!/' 3/2 DTF yields consistently higher FFs. 


very close to the Newtonian expression As a result, the improvement [(■i/'o ^ 3 / 2)2 to (^/>o '03/2 0 ^) 4 ] in FF is 

only ~ 0.3-1.6%, while (at least according to the simple Gaussian analysis of Sec. IV Bll the detection threshold 
increases by ~ 4% (although this number does not take into account the 0“ = 0 constraint). It seems therefore 
that the (00 03/2 a )4 DTF is not a useful upgrade of (0o 03 / 2)2 for the purpose of detecting the signals emitted 
by precessing binaries. 

4. The (00 03 / 2 ' 6)0 DTF [lighter-red solid lines] includes modulational corrections for both amplitude and phase. 

The resulting improvement in FF over the SPA families is remarkable (for BBHs, 8-22% over SPAs, and 6-10% 
over SPAc; for NS-BH binaries, 20% over both). However, the effect of the modulational terms is seen best by 
comparing (0o 03/2 B)q to (0o we get an improvement of 5-9% for BBHs, and 20% for NS-BH binaries. 

This numbers should be compared with the projected increase ~ 8% in the detection threshold (Sec. |vU. 

5. For the (0o 03/2 DTF, where the frequency dependence of the modulating terms is f~^l^ rather than f~^, 
fitting factors are not significantly different from (0o 03/2 B)e. Therefore we do not show these numbers. Tables 

Enn and im also contain a few FFs computed against the 3.5PN order target model (with 9 = 0). The FFs, 
shown for the (0003/2^)6 DTF, are essentially in line with their 2PN counterparts. 

Our results suggest two strategies to search for the signals from the precessing BBHs examined in this paper. We 
can try to follow the modulations induced by precession, using a DTF similar to (00 03/2^8)6; or we can just use 
(00 03/2)2, which is considerably better than SPAs (mostly because of /cut), and slightly better than SPAc (because 
of the extended parameter range). The gain in FF when we upgrade (0o 03 / 2)2 to (00 03/2 8)6 is offset by a similar 
increase in the detection threshold, but the latter increase might be contained by reducing the range of the allowed 
ctfc, or by other data-analysis considerations that do come into the simple Gaussian analysis of Sec. IvH 

Figure IT^ shows the projection of the 2PN target waveforms onto the (0o, 03 / 2 ) section of the (0o03/2)2 parameter 
space; Fig. ca shows the projections of the waveforms onto the ( 0 o, 03 / 2 ) and (0o,8) sections of the (00 03/2 8)6 
parameter space. It is interesting to notice that, with either strategy, the ranges of 0o and 03/2 needed to match 
effectually the signals in our populations are essentially the same found in BCVl to match the signals predicted by 
a variety of PN models for BBHs without spins. In Figs. d and d these ranges are delimited by the thick dashed 
lines; the thin mass lines represent the range of detection templates needed to match effectnally the signals predicted 
by different PN models for the same binary masses. As we can see, the projections of the spinning-binary signals are 
smeared around the nonspinning-binary mass lines with the same masses. 
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FIG. 12: Projection of the 2PN target signals onto the {'ipo'4’3/2)2 DTF. For the (10 + 10)Mq, (15 + 15)Mq, (20 + 5)Mq, 
(7 + 5)Mq, and (20 + 10)Mq BBHs in our Monte Carlo populations, the clusters of gray dots show the projection of the 2PN 
target waveforms onto the (i/>o,V’ 3 / 2 ) parameter plane of the {tpo4’3/2)2 DTF (the projection of a given target signal is given 
by the values of tpo and i/’ 3/2 that maximize the FF; here /cut is not shown). For each set of masses, we draw a dashed ellipse 
centered on the parameter-space baricenter of the dots, and sized to include 90% of the dots (the proportions of the axes follow 
the two-dimensional quadratic moments of the dots). The larger blue dots, joined by the thin lines {mass lines), show the 
projections of the nonspinning PN models studied in BCVl, for the same sets of masses plus (5 -I- 5)Mq and (10-1- SjM©; each 
line joins signals with the same binary masses, but obtained from different PN target models. As we can see, for each set 
of masses, the projections of the spinning-binary signals are clustered around the corresponding mass line; moreover, all the 
projections fall within the region (delimited by the thick dashed lines) suggested in BCVl to match all the nonspinning PN 
models. 


Thus, a signal search based on the {ipo'^ 3 / 2)2 DTF is a good starting point for both nonspinning and spinning 
binaries. It might also pay off, depending on the results of a more realistic evaluation of false-alarm probabilities, to 
upgrade this DTF to (V’o V^ 3/2 ct)4j with improved performance for nonspinning but nonadiabatic BBHs, as shown in 
BCVl; or even to (V’o ' 4 ^ 3/2 B)q, with the best FFs for spinning binaries and without any deterioration for nonspinning 
ones. 


D. Modulated DTFs for NS—BH binaries 

Let us now look in detail at the FFs achieved by the DTFs and standard template families against the signals 
generated by (10-1- 1.4)M0 NS-BH binaries where the BH is spinning rapidly (see Tab. lIXI and Figs.lllland lldll . First 
of all, we notice that there is little difference between the performance of the SPAs and SPAc templates, because the 
ending frequency lies outside of the band of good interferometer sensitivity. Furthermore, the number of GW cycles 
within this band is very high, so it is crucial that a DTF reproduce very accurately the evolution of the GW phase; so 
using the {tpo 143 / 2)2 DTF improves only slightly on the performance of the SPA templates. Introducing precessional 
corrections brings about a dramatic change: for the {ipo 143/2 DTF, the increase in FF and FFes with respect to 
SPA is respectively 20% and 16%, which is enough to justify the introduction of six ak coefficients, according the 
Gaussian analysis of Sec. EH 

The dependence of the FF on the spin configuration is shown in Fig. ^1 For the NS-BH signals in our Monte 
Garlo population. Figs. ESI and d show the template parameters i 4 o, 143/2 j and B that maximize the overlap plotted 
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FIG. 13: Projection of the 2PN target signals onto the (V’o '03/2 Sje DTF. For the (10 + 1O)M0, (15 + 15)Mq, (20 + 5)Mq, 
(7 + 5)Mq, and (20 + 1O)M0 BBHs, and for the (10 + 1.4)M0 NS-BH binaries in our Monte Carlo populations, the clusters of 
gray dots show the projection of the 2PN target waveforms onto the {ipo, 03/2) [on the left] and (0o, B) [on the right] parameter 
plane of the (0o03/2 B)^ DTF. For each set of masses, we draw a dashed ellipse centered on the parameter-space baricenter of 
the dots, and sized to include 90% of the dots (the proportions of the axes follow the two-dimensional quadratic moments of 
the dots). The nonspinning-model mass lines and the boundary of the suggested parameter ranges are shown as in Fig. 1121 



/^ax 
''■eff ''■eff 


FIG. 14: Average fitting factor for the DTFs and for the SPAs template families for (10 -I- 1.4)M0 NS-BH binaries, plotted 
against the initial Kefr = Ljv • Seff. The vertices of the segmented curves show the FF averaged on the sets of samples that 
fall within the equispaced Kefr bins [—1, —0.8), [—0.8, —0.6), ..., [0.8,1]. The error bars show the sampling error on the bin 
averages. We plot also two additional vertices, aligned with the abscissae —1 and 1, which show the FF averaged over the 
bins [-1,-0.98) and [0.98,1]. 


against the initial k (conserved in NS-BH binaries). In the left panel, we see that the parameter 0O) which is related 
to the Newtonian chirp mass, has only a weak dependence on k (it varies by ~ 8%); on the other hand, the parameter 
03/2 has a strong dependence. A plausible explanation is that the SO term in the SPA phasing is formally 1.5PN 
[see Eqs. and ra i. and so is the term 03 / 2 / in 0 nm(/), which takes on the job, as it were, of reproducing the 
nonmodulational effects of the SO coupling. In the right panel, we see that for most of the binary configurations the 
values of B cluster around three lines [B = 100, B — {\ + n) 110 -I- 110, and B — {1 + k) 240 + 160]. Further analysis 
are needed to provide an explanation for this interesting behavior. 

Thus, the (0o 03/2 B)e DTF is a good candidate for the data-analysis problem of detecting GW signals from NS-BH 
binaries with rapidly spinning BHs. However, the analysis of processional dynamics and GW emission carried out in 
this paper suggests an even more specialized DTF, which could be built with the following guidelines. 

1. The waveform can be computed directly from Fq. (18211 (obtained in the processing convention): the necessary 
ingredients are the time evolution of the orbital phase fk and of the binary polarization tensors [e+^x]i/, plus 
the fixed detector polarization tensors [T+_x]i/- 
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FIG. 15: Projection of the (10 -I- 1A)Mq NS-BH target signals (computed at 2PN order) onto the (^/jq 'i/’ 3/2 S)6 DTF. The dots 
show the values of the V'o (left panel) and 4 ’ 3/2 (right panel) target parameters that yield maximum overlaps with the signals 
in the target populations. 


2. The evolution of ili is obtained by solving Eq. o, where S 2 can be set to zero, and Si enters only in the 
conserved term L^v • Si. As a consequence, Eq. is effectively uncoupled from the evolution of L^v, Eq. 

3. The evolution of the tensors [e+_x(t)]y is obtained from Eq. 17111 . after integrating Eqs. and for the 
coupled evolution of La^ and S, which depends only on Lat • Si, on Si (conserved), and on w(t). 

4. A source frame attached to the initial configuration of the binary, similar to the frame constructed in Sec. IIV Al 
[see Eqs. mi can be used to carry out the explicit construction. By way of the initial conditions 

the tensors e_|_ ^ and the orbital phase di (up to an additive constant Tq) are then well defined as functions of 
the basic and local binary parameters only. We have therefore completed the specification of the first part of 
Eq. 1X211 . which expresses the components of the mass quadrupole moment. 

5. The remaining part of Eq. (EH), which expresses the projection on the polarization tensor of the detector, 

pv = [T+rF+ + [T^Y^F^, (120) 

is determined by the directional parameters 0, ip, (j), 0, and ip, which are now referred to the source frame 
attached to the binary. When we look for GWs using matched filtering, we can search rapidly over such a 
parametrization by treating the P*-' as extrinsic parameters, along with the time of arrival and the initial orbital 
phase The only intrinsic parameters would then be toi, m 2 , Si, and S • La/, all of which are conserved. 

This family of templates adds a further intrinsic parameter with respect to {ipo'4’3/2^)6, but it has the advantage of 
producing essentially exact waveforms (valid in the adiabatic regime, and up to the highest PN order included), and 
of expressing these waveforms directly in terms of the physical spin parameters Si and S • La/. We believe that the 
implementation and the false-alarm statistics of this family are worthy of further investigation |39|. 


VII. SUMMARY 

In BCVl, the nonmodulated DTFs (V'oV' 3 / 2)2 and {ipQip 3 / 20 t)i were shown to have FF > 0.95 against several 
nonspinning-BBH target models, obtained under different PN approximation schemes. In this paper, we have shown 
that these two families are also rather effectual at matching the signals from BH-BH and NS-BH processing binaries 
with single-BH masses between 5 and 20 Mq and with maximal BH spins, at least if these signals can be described 
by an adiabatic sequence of quasicircular orbits up to 2PN order. 

More specifically, for (7 -I- 5)Mq, (10 -I- 1O)M0, (20 -I- 1O)M0, and (15 -I- 15)M0 BBHs, we obtain FF > 0.93 
and FFeff ^ 0.95. The improvement is 2-16% over Schwarzschild-terminated SPAs templates, thanks largely to the 
ending-frequency parameter /cut; and 1-2% over SPAc templates, thanks to the effective extension in the range of 
parameters, released from their functional dependence on the masses of the binary. Although the latter improvement 
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FIG. 16: Projection of the (10 + 1.4)M0 NS-BH target signals (computed at 2PN order) onto the {ipo i/’ 3/2 B)e DTF. The dots 
show the values of the B target parameter that yield maximum overlaps with the signals in the target populations. 


seems negligible, we should keep in mind that ('0o'03/2)2 DTFs are also more suitable to match the nonspinning 
BH binaries studied in BCVl with PN expanded and resummed models. Results are worse for binaries that have 
smaller mass ratios rj, and therefore more GW cycles in the band of good interferometer sensitivity. In this case the 
modulational effects due to precession become important, and must be included in the detection templates. Indeed, 
for (20 + 5)Mq BBHs, the { 1 ^ 01 ^ 3 / 2)2 and (0o'03/2a)4 DTFs have FF ~ 0.89 and FFeff ~ 0.92; for a (10 + 1.4)M0 
NS-BH binary, we find FF ~ 0.78, and FFes ~ 0.83. 

Motivated by these shortcomings, we have investigated in detail the dynamics of precession in these binaries, and 
we have introduced a new convention to write the GW signal (as computed in the quadrupole approximation) as 
a function of binary and detector parameters, isolating the oscillatory effects of precession in the evolution of the 
polarization tensors [e+_x]i/- As a result, the detector response to GWs can be written as the product of a carrier 
signal, which very closely resembles the nonspinning signals studied in BGVl, and a modulational correction, which 
can be handled using an extension of Apostolatos’ ansatz (IMI) . On the basis of these observations, we build the 
modulated DTF (00)'03/2^8)6, which yields FF and FFeff — 0.98-0.99 for the BBHs investigated, and FF ~ 0.93, 
FFeff — 0.95 for (10 + 1.4)M0 NS-BH binaries. This DTF has the advantage that all the modulational parameters 
(except for B) can be treated as extrinsic parameters, reducing considerably the computational cost of signal searches. 
According to the simple analysis of Sec. IvH the detection thresholds for this DTF should be set higher than those 
for simpler families; still, the gain in FF is still somewhat larger than the increase in the threshold, and more realistic 
analyses of false-alarm statistics might provide a way to sidestep this difficulty. The same arguments that lead to the 
(00) 03/2'S )6 DTF suggest a new, very promising class of templates for NS-BH binaries, which we discuss briefly in 
Sec. IVI l51 and which we plan to investigate more thoroughly elsewhere 

We wish to make a few final remarks. First, in this paper we limited our analysis to compact objects moving on 
quasicircular orbits; from the results on the ending frequencies (see Fig. 0 we see that there exist spin initial conditions 
for which the ending frequencies (end of inspiral) are in the LIGO-VIRGO band. So, in these cases we should use 
spinning dynamics that goes beyond the adiabatic approximation. This dynamics (without radiation-reaction effects) 
is already available in the FOB framework ill thanks to the work of Damour M- We plan to investigate the effects 
of nonadiabatic PN dynamics in the near future. 

Second, a few years ago Levin pointed out m that spin-spin effects can introduce chaos into the trajectories; as 
a consequence, the gravitational waveforms would come to depend sensitively on the initial conditions. More studies 
followed Considering only conservative dynamics (no RR), Cornish and Levin m found some examples of 

rather eccentric (e ^ 0.6 or 0.9) chaotic orbits, and a few quasicircular chaotic orbits. However, these authors observed 
that chaos would be damped by RR effects, and that it would not affect the inspiral waveforms, except (perhaps) 
at the very end (the plunge). Still, at this time the dynamical structure of phase space has not been explored 
systematically, and a more conclusive study tuned to the LIGO-VIRGO detection problem remains desirable. The 
analysis of this paper assumes that, by the time the GW signal enters the band of good detector sensitivity, RR effects 
have circularized the orbit, and have brought the binary into the adiabatic regime, which is valid until the MEGO. We 
did not try to perturb the initial conditions slightly and to investigate the resulting changes in the orbital evolution 
and in the waveforms. 

Third, we have evaluated the performance of our DTFs by averaging over uniform distributions of the initial spin 
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angles. Of course it would be preferable to assume more realistic, nonuniform distributions derived from astrophysical 
considerations. Some results for spin distributions in BBHs (with only one spinning BH), and in NS-BH binaries 
were obtained by Kalogera using population-synthesis techniques ji^. In particular, Kalogera found that 30-80% of 
the NS-BH binaries that will coalesce within a Hubble time can have a tilt angle (the angle between the spin and the 
orbital angular momentum) larger than 30 degrees. These results assume that the spinning BH in the binary forms 
first, and that its spin is aligned with the orbital angular momentum; the tilt angle originates from the supernova 
explosion that forms the NS. For the case of the binaries formed in globular clusters, there is no theoretical argument 
to suggest any particular spin distribution. 

Finally, recent analyses of spin-spin effects in the PN inspiral equations M suggest that, for comparable-mass 
BBHs, by the time the GW signal enters the band of good interferometer sensitivity the two BH spins may have 
become roughly locked into a fixed relative conhguration. If these results are confirmed, they could provide preferred 
initial spin conditions, and simplify the data-analysis problem for comparable mass binaries, by reducing the variability 
of expected GW signals. 
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APPENDIX A: VALIDITY OF THE ADIABATIC SEQUENCE OF SPHERICAL ORBITS 

In the target model defined in Sec. inn the inspiral of the two compact bodies is described as an adiabatic 
sequence of spherical orbits. In this Appendix we wish to discuss the validity of this assumption. Introducing the 
orthonormal basis (A, n, Ljv)) where n = x/r, L^v = Ljv/Tat, A = L^v x fi and L^v = /ix x v (with ^ the reduced 
mass), it is straightforward to write the equations of motion as [see Eqs. (4.1) of Ref. we use the relations 
V = r n -I- r w A, r'^ 


n • a = 
A • a = 

Lat a = 

where a is the acceleration in harmonic gauge given b 
Eq. ESI) then implies w = 0; and from Eq. m we get 


r — ruj'^ , 

(Al) 

ruj + 2ruj, 

(A2) 

dliN ~ 
dt 

(A3) 


Eqs. (2.2a), (2.2c) of Ref. [T^. If we impose r = 0 = r, 


2 2 ^ 
r UJ = - 
r 



So; = 


2 mi / 


V 2 m2. 


(A4) 


where for simplicity we have set M = 1. Although spherical orbits (orbits where both r and to remain constant) exist 
at any given instant, they are not preserved along dynamical evolution because the quantity Ljv • that appears in 
Eq. 1A4II is not conserved. Indeed, averaging over an orbit (and, for simplicity, neglecting spin-spin effects), we 
get 


dt / 


(A5) 
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where (A) denotes the quantity A when the spin-orbit (and spin-spin) terms have been averaged over an orbit. Using 
the precession equations for the spins we derive 


d(L 


N 


dt 


= —6 - Ln 

mi m 2 


(Si X 82)4.. 


(A 6 ) 


Hence, because the circular-orbit condition is not preserved during the evolution, either (w) ^ 0 or {r)^ 0 (or both). 

Let us now see how Eq. o for oj changes if effects of this kind are included. The usual argument [la [^ used to 
obtain the adiabatic evolution of uj rests on the energy-balance equation. 


Arr = —E{uj, Lat, Si, 82 ) = 


dE . f dE 


dE ■ dE 


duj'~ ' VSLat ^^”^981 ^^”^582 * 


(A7) 


where 


A(a;,Liv,Si, 82 ) = -|(Ma;)2/3| ^ 


3M2 


Lat • 8 eff (Muj) ■ 


(AS) 


^(_81 + 57 rj + ^ [(81.82) - 3(LAr • Si)(Liv • S2) 




is the orbital energy evaluated at Newtonian order, but including spin-orbit and spin-spin effects, and where Arr 
is the RR energy loss laiii. From Eqs. (TASll . ®, @ and 0 , we notice that the sum of the last three terms in 
parentheses in Eq. does not vanish: at leading order, its value is 


-£/fixtra, — 




81 X 82 -L 


-‘N 


(A9) 


This expression is zero if masses are equal, or if spins are either aligned or antialigned. Retaining the term (IA9II in 
the calculation yields an additional contribution in the evolution of cj, with a leading order correction 


(^extra 



3 (toi — m 2 ) 

4 M 


VXiX 2 iMuj)'^ (^81x82) 



(AlO) 


Thus, compared with the other terms in Eq. o, Wextra appears formally at 0.5 PN order (very low!) in the expansion 
of w. Note that the spin-orbit term in the energy (IA 8 II . combined with the leading-order precessions, does not produce 
such a term; this makes the adiabatic approach fully consistent up to 1.5PN order. In fact, Eg^tra originates from 
taking the derivative of Ego and using next-to-leading order terms in the precession equations, and the derivative 
Egg while using the leading-order terms in the precession equations. 

However, the effect of this term in the regime that we consider is not as large as suggested by its formal PN order. 
For example, under the worst possible assumption (that the geometric factor [(81 x 82 ) -Lat] has always the maximum 
value of one, and that spins are maximal), we get the correction 


extra 

2tt 


1 25 

2t: 16384 rj 


(Ma;/)-4/3 




(All) 


to the number of orbital cycles, where LUi and ujf are the initial and final orbital frequencies under consideration. This 
is formally a 0.5PN correction, as can be seen by comparing it with Eq. (4.16) of Ref. ^3. Nevertheless, for (say) a 
(20 -I- 5)Mq BBH, this correction will be at most 0.34 orbital cycles from Wi = tt x 30Hz to Wi = tt x 400 Hz, to be 
compared with a baseline of 52 orbital cycles from the Newtonian term and 8 from the IPN term. For a (10-I- 1 . 4 )M 0 
binary, the correction will be 1.6 orbital cycles, to be compared with 175 orbital cycles from the Newtonian term 
and 30 from the IPN term. The correction is small because, although the PN order is formally low, the numerical 
coefficient of the geometric factor [(Si x S 2 ) • LAf] is very small. 

So far, we have assumed [(Si x 82 ) - La/] ~ 1 along the evolution. Let us now estimate the more important effect 
that comes from the precession of Lat, Si and 82 , which is especially important for binaries with small mass ratios, 
which have longer RR time scales and more precessional cycles. At the leading order (with M = 1) 


_d 

dt 




(mi - m 2 ) [Si • 82 - (Si • Liv) (S 2 • Liv)] + 0 ( 0 ;^), 


(A12) 
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and 


£ 

dt 


Si-S2-{Si-Ln){S2-Ln) 


-{mi — m2) ^81x82^ 


L 


N 


+ 0{u;^). 


Combining the above equations, we get (at leading order) 


(A13) 


dt"^ . 


81 X 82 


) 






(A14) 


This means that the geometric factor [(Si x 82 ) • Lat] oscillates around zero with a time scale ~ . Thus the effect 

of Wextra accumulates only within this timescale, which is 1.5 PN orders shorter than the RR timescale. Therefore, 
we expect that the real extra will be even smaller than the formal prediction given by Eq. mn), and that it will 
contribute effectively at 2PN order. As a check, we evaluated the FF between the gravitational waveforms obtained, 
for a (10 + 1.4)Mq binary, by first including and then dropping the extra term in Cj. We found that the FF is ~ 0.99. 
On the basis of this last check and of the analysis outlined above, we conclude that the adiabatic assumption is quite 
adequate for the purposes of this paper. 


APPENDIX B: PROOF THAT THE PRECESSING CONVENTION YIELDS u = i>s 

First of all, it is easy to confirm that, as long as 61^2(0) and Ljv(O) form an orthonormal basis at some initial time, 
the evolution equation ei^2 = x ei_2 will always keep the triplet an orthonormal basis. It is then always possible 
to have a 4)(t), such that 


h(t) = ei cos$(t) + e2sin$(t), A(t) = —ei sin<i)(t) + 62 cos$(t). (Bl) 

Taking the time derivative of h(t), we have 

h = $A + Oe X n, (B2) 

Now, the adiabatic condition for a sequence of circular orbits states that h = wA, so we have 

n = wA = $A + Oe X h. (B3) 

By definition [Eq. iTT^ j. fig has no components along 63 = Lat. It also has no components along A, because 

Jig X Liv = Lat = nxA + nxA = a;AxA + hx (—<i)n + Jig x A) (B4) 

= r2e(n • A) — A(n • Jig) oc A, (B5) 


where in the last step we used h • A = 0 and the vector-triple-product rule. It follows that Og lies along n, and 
therefore Hg x n = 0. Equation then gives the desired result, = w, i.e. $(t) = vl/(t) -|- const. 
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